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SUMMAEY 


An accurate method of determining secondary stresses in thin- 
walled, uniform beams of closed cross section is herein presented. 

The cross sections are assigned to he preserved hy closely spaced rigid 
diaphragms. In section I the integrodifferential equation governing 
axial displacements is formulated and solved for a beam without longi- 
tudinal stiffeners. In section II the corresponding summation- 
difference equation is developed and solved for a beam with stiffeners 
(flanges and stringers) . The cross section, loading distribution, and 
end conditions are assinaed to be arbitrary. 

By introducing generalized difference equations the mathematical 
analysis for the stiffened beam may be performed in a manner exactly 
analogo;is to the process used for the unstiffened beam. A separation 
of variables in the homogeneous eqiiation leads to the natural stress or 
displacement modes for a cross section. The solution of the nonhomo- 
geneotis equation is then expressed as an expansion in terms of the 
natural stress modes. Particular attention is given to cross sections 
with single symmetry and double symmetry. 


HISTOEICAL NOTE 


Eeference will be made only to papers which give exact solutions 
for beams having rigid bulkheads with monocoque ( unstlf f ened) or semi- 
monocoque (stiffened) sections. In most cases the solutions for s emi - 
monocoque sections are those obtained by replacing the true section by 
an Idealized section in which the normal stresses are carried on a 
finite number of stiffeners and uaiformly distributed shear flows are 
carried on the connecting webs or wall elements. 

If an idealized cross section has only three stiffeners, or flanges, 
the distribution of axial stress must be planar and there are no sec- 
ondary stresses. The doubly symmetrical rectangular idealized section 
with fo\ir flanges was intreduced and analyzed by Ebner in 1933 
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(reference l) . This cross section lias one secondary stress mode "wlilcli 
occurs under torsional loading. Ebner considered such, modifications of 
the problem as the effects of treh buckling, shear deformation of the 
ribs, and finite rib spacing. The solution for the four- flange doubly 
symmetrical section was also determined independently by Grzedzielski 
in 193 ^ (reference 2) This author recognized the Importance of 
Idealizing the cross section for stress-analysis purposes. A con^iari- 
son of the solution for the four- flange section in torsion with eitperl- 
mental tests on a box beam was given by Kuhn in 1935 (reference 4). 

A four-flange section with curved webs of arbitrary shape was 
analyzed by KLrste in 1937 (reference 5) . This author gave the loca- 
tion of the principal shear axes for any arbitrary cross section. He 
also gave the formulas for the associated shearing section properties. 
The torsional section constant was called the central moment of inertia, 
while the areas resistant to horizontal and vertical shear were called 
the reduced areas. This author also located the zero-warping axis for 
a four- flange beam. He showed that waorping displacements in a four- 
flange beam are due entirely to the torque about the zero-warping axis. 

The first exact solution for a secondary stress mode in bending 
was given by Kuhn in 1937 (reference 6). The cross section considered 
by him was a rectangular, five-flange, open section in which the bot- 
tom wall was missing. The section was symmetrical about a vertical axis 
and had a single central stringer on the top wall. Experimental test 
results were compared with the theoretical solution. 

A general solution for stiffened beams of arbitraxy cross section 
was given by Ebner and Koller in 1937 (reference 7) • The beam was 
assumed to have rigid or elastic stiffening rings at finite spacing and 
to carry bending or torsional loads. The secondary self-equilibrating 
stress distributions were chosen in an arbitrary manner. Energy prin- 
ciples were iised to determine the spanwise variation of the secondary 
stress distributions. The underlying viewpoint of the analysis con- 
sisted of regarding the structure as a space framework. The calcula- 
tion effort required by the method is so extensive as to make the method 
iB 5 )racticable except for a beam having a small number of stiffeners and 
stiffening rings. The specific cross sections considered in detail by 
Ebner and Kol-ler were elliptical sections with fo\ar stiffeners and with 
six stiffeners having double symmetry. They also considered a circular 
section with twelve stiffeners and double symmetry. These authors did 
not introduce the principal shear axes or the associated section 


■^Jhe 


formulas of reference 2 are derived in a different manner in 
reference 3 where a corrected formula is given for the rate of twist at 
the root of a cantilever beam. 
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properties. A comparison was given "between tlieo27y and experiment for 
a beam of circiilar cross section. 

Following the method of analysis that had been given by Kirste, a 
more comprehensive treatment of the four-flange beam of arbitrary cross 
section was given by Drymael in 19^1 (reference 8) . This author gave 
many new formulas and equations for the three-flange section and the 
four-flange section. He gave a careful, detailed study of the geometric 
properties of the shear-carrying area of idealized stiffened sections. 

He showed that the "reduced areas" are a maximum and a minimum with 
respect to the principal shear axes. He also showed that the "central 
moment of inertia" is a m1 n iT rnim when computed with respect to the prin- 
cipal shear axes. Formulas applicable to arbitrary cross sections were 
given for lateral deflections due to shearing strains in thin-walled 
beams. Strangely enough, these formulas, even today, are not well- 
known even though the deflections due to shearing strains are of appre- 
ciable lB5>ortance in aeroelastlc and dynamic analysis of airplanes. 

TTie first exact analysis for a beam of monocoque, or unstiffened, 
cross section was given by "Von Karman and Chlen in 19^ (reference 9) • 
Their solution is limited in its applicability to sections having double 
symmetry. They considered only fhe torsion case and assumed the beam 
to have closely spaced rigid diaphragms. The integrodifferential equa- 
tions governing the axial displacements were formulated and solved for 
a semi-infinite cantilever beam acted upon by a torque at the tip. 

Stress distributions were . obtained for a rectangular section and a 
rectanglelik;e section with rounded comers. 

A very general and comprehensive exact analysis of beams with 
rigid diaphragms was given by Hadji-Argyris and Dunne in 19^7 (refer- 
ence 10) . The end conditions, loading conditions, and cross section 
were considered to be arbitrary. Conical as well as \miform beams were 
analyzed with multicell or single-cell cross sections. Both open and 
closed cross sections were considered with and without stiffeners. It 
is impossible to summarize the contents of the paper here because of 
its great length. The method of analysis used by these authors is 
entirely different from that employed by "Von Karman and Chlen. The 
differences between the methods are discussed in the text. These 
authors found a zero-warping axis to be associated with each secondary 
stress mode corresponding to the axis found by'Kirste for the four- 
flange section. They illustrated clearly the large difference which 
may exist in the positions of the zero-warping axis and the shear center 
for a four-flange beam. A detailed analysis of a four- flange tapered 
wing with cut-outs was presented. 

Specific formulas were given by Kempner in 19^7 (reference ll) for 
a six-flange doubly symmetrical hexagonal section. A five-flange open 
section, obtained by removing the bottom webs and bottom stringer, was 
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also considered. The beam was subjected to vertical bending loads only. 
Theoretical solutions were congjared with experimental results which had 
been previously published by Peterson in 19^5 (reference 12). 

A rather sin^jle method of analysis was given by Levy in 19^7 
(reference 13 ) for beams with arbitrary stiffened cross sections but 
limited to beams with finite bTilthead spacing. The loading was assmed 
to consist of concentrated loads applied at the intersections of the 
ribs and spars. The solution involved the use of eguilibrium conditions 
and an energy principle. The method is applicable to beams with cut- 
outs and sweepbach. However, a comparison of theoretical solutions by 
Levy’s method with experimental results, for a swept beam, was given by 
Bisplinghoff in 19^9 (reference l4) and did not show satisfactory agree- 
ment. Deflections rather than stresses were compared. In the theory of 
Levy and the beams tested by Bisplinghoff the ribs of the swept beam 
were parallel to the plane of s ymme try of the airplane. 

The solution for a swept beam with closely spaced diaphragms was 
given by Wittrick in 19^ (reference 15) • The method of analysis used 
was an extension of the theory of Hadji-Argyris and Dunne (reference lO) 
to Include the effect of sweepback. The general integrodifferentlal 
equation for the stress function for a beam of arbitrary cross section 
was formulated. Specific solutions were given for uniform and 
exponentially varying section torque and bending moment. 

An exact solution for a single-cell monocoque (unstiffened) beam 
with closely spaced diaphragms was given by Adadurov in 19^ (refer- 
ence 16 ) . The cross section and loading were assumed to be arbitrary. 
This author introduced a stress function having its first derivative 
proportional to the axial displacements and second derivatives propor- 
tional to the stresses. The method is actually a minor modification of 
the method of Von Karman and Chien (reference 9) but is extended to 
arbitrary cross sections and loading. The author’s method of solution 
is unwieldy because he failed to introduce the principal shear axes. No 
specific solutious were obtained. 

Two specific solutions for swept beams were given by Thompson and 
Wittrick in 19^9 (reference 1?) using the theory that had been pre- 
viously published by Wittrick (reference 15) • The first solution was 
given for a singly symmetrical trapezoidal section with four flanges 
and closely spaced diaphragms located in a skewed position with respect 
to the axis of the beam. The walls were assumed to carry both normal 
and shearing stresses. The loading consisted of a concentrated moment 
and a concentrated torque at either end (the St. Venant problem). In 
the second example the cross section was a doubly symmetrical rectangle 
and the loading was assimied to produce a section torque varying 
exponentially. 
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A general analysis of stiffened Toeams witti artitraiy cross sec- 
tiona vas given ty Duberg in a thesis siibmitted to the University of 
Illinois (Feh. 19^). An abstract of this thesis vas published sub- 
sequently by Duberg in 19^9 (reference I8) . The problem treated by 
this author may be regarded as an extension of the four-flange problem 
treated by KLrste and Drymael to sections with more than four stiff- 
eners. The method of solution used by Kirste and Drymel cannot be 
extended directly to more than four flanges and new methods of solution 
are reqiiired. Duberg 'a method of solution resembles that of Von K^mian 
and Chien in that the axial displacements were chosen as the unknown to 
be determined. Princix)al shear axes play an essential role in this 
method of analysis and convenient formulas were given for their calcu- 
lation. The formulas correspond closely to those given by Kirste and 
Drymael. numerical methods were employed to obtain two solutions for 
cantilever beams having rectangular cross sections with single symmetry. 
Both finite 2rib spacing and closely spaced ribs were considered. 

The method of analysis to be given in this paper may be considered 
as an alternative to the method of Hadji-Argyris and Dunne. It corre- 
sponds to, and employs ideas contained in, the papers by Kirste, 

Drymael, Von Karmi^ and Chien, and Duberg. 

This work was originally a thesis for the degree of doctor of 
philosophy at the California Institute of Technology and has been made 
available to the MCA for publication because of its general interest. 


I - BEAMS HAVING THIN-WALLED UNSTHFENED CEOSS SECTIONS 

Intr oduc t ion 


In section I the beam will be assumed to have a thin-walled, 
■unstiffened cross section. It -will be assumed to consist of a single 
cell -without comers. This is the simplest and most convenient section 
to consider in developing a rational theory for the determination of 
secondary stresses in thln-'walled beams of arbitrary cross section. The 
loading conditions and end conditions for the beam nay be of any 
arbitrary nature. The beam is assumed to be of uniform section -with 
no cut-outs. 

The -wall of the cross section is assumed to be sufficiently thin 
that one imy consider the stresses to be uniformly distributed o"ver the 
thickness of the •wall. The thin shell then acts essentially as a 
membrane. The shell is assrmied to be stabilized by closely spaced 
rigid diaphragms -which preserve the cross-sectional shape of the beam 
under the action of loads. It -will be assumed throughout the analysis 
that the diaphragms are infinitely stiff in their own planes but 
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completely free to warp out of their planes. This asstmqotion is the 
only one ■which has a significant effect ■upon the final results and is^ 
th^us, the principal source of any error which may he contained in the 
solution. 

Upon replacing the real, physical heam hy a similar heam ha^ving 
rigid diaphragms, one obtains a ■well-defined problem in stress analysis 
for which exact soltrtions can he oh"bained. Such exact solutions of the 
similar heam may, in most cases, he regarded as "acc^ura'te" (hut not 
exact) solutions for the real, physical heam. One may also obtain 
approximate solutions for the similar heam "with rigid diaphragms. Such 
solutions may he regarded as "approximate" solutions for the real 
physical heam. l&ny such approximate solutions ha^ro been piiblished for 
beams ha'vi.ng sln5)le cross sections. 


Symbols 


A 

Aw 

% 

Ay 


E 


cell area 

area of monocogue cross section 
horizon'tal shear-re sls'bant area 
vertical shear-resis^tant area 

section cons^tants associated -with ith stress mode 
Young's modul'us 



f 

F 

g 

G 

h 


transverse distrih\rtion of axial displacements or normal 
stresses 

stress function 

span^wise -variation of axial displacements 
shearing modul-us of elasticity 
spanwlse variation of norml stress 


H 


horizon'tal section shear 
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Ic 

Ix 

J 

Mx 

“y 

p 

% 

Pv 

p 

q 

r 

s 

t 

T 

u 

V 

V 
w 


x,y,z 


x,y 

a 


central moment of inertia 

polar moment of inertia 

moment of inertia about x-axis 

moment of inertia about y-axis 

torsion constant 

bending moment about x-axis 

bending moment about y-axis 

normal-stress flow in monocoque sections 

applied horizontal loading 

applied vertical loading 

axial force in beam 

shear flow 

radius to a tangent 

tangential coordinate 

applied torsional loading 

wall thickness of monocogue sections 

section torque 

horizontal displacement 

vertical displacement 

vertical section shear 

axial displacement 

rectangular coordinates 

centroidal coordinates 

angle between tangent and x-axis 
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e 

0 


X. 


= 



sliearirig strain 
normal strain 

angle "between printipal bending axes and principal shear 
axes 

ch£iractei*istlc number 
also Poisson's ratio 


cp angle of twist (rotation) 

I tangential displacement 

CT normal stress 

T shearing stress 


Fundamental Equations of Elasticity 

The beam and coordinate system are shown in figure 1. Boundary 
conditions on stresses or displacements must be Tmown at each end of 
the beam. A segment of the curved wall is shown in figxire 2(a) . It is 
convenient to introduce the coordinate s which is measured aromd the 
periphery of the beam along the center line of the wall. A point in 
the wall of the beam may be located by givirg the values of s and z 
rather than the values of x, y, and z. The coordinate s and the 
shear flow q are considered to be positive in the coimterclockwlse 
direction as shown in figure 2(a). 

The differential element dsdz is shown in figure 2(b) with the 
forces which act on it in the axial direction. The remaining forces 
which act on the element are shown in figures 2(c) and 2(d). Assuming 
that Hooke’s law is applicable, the stresses and strains are related by 
the following well-known equations: 

- ^^s = (1) 


- ^z= ®^s (2) 


^sz ~ ^^sz 


( 3 ) 
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Because of the assTm^jtion of infinitely stiff diaphragms^ one may 
immediately set eg equal to zero, since such diaphragms would prevent 
any strain in the wall in the tangential direction. Using equation (2) 
the following equation is obtained: 

= 0 

Og = p0z 



From this equation the tangential normal stress may he computed at any 
point after the axial normal stress distribution has been determined. 
Equation (4) may be substituted in equation (l) to obtain cfz in terms 
of gz* 


a 


z 



(5a) 


where 


E* 



(5b) 


The strains at a point may be expressed in terms of the axial dis- 
placement w and the tangential displacement | by the following 
well-known equations: 


^z = 




(6a) 


7 


sz 


^ hi 


(6b) 


The axial displacement is measured as positive in the positive direction 
of the z-axis while the tangential displacement is positive in the posi- 
tive direction of the s coordinate, that is, counterclockwise. The 
stress Ss a.nd the strain 69 do not appear in the analysis henceforth 
and the subscripts of the remaining stresses and strains will be dropped 
with the following change of notation: ^ ~ ® ~ ^ ^sz^ 

/O. - _ 
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elasticity gives the following formulas for stresses in terms of 
displacements: 


0- = E‘ 




(7a) 


T 



+ G 


hi 


(7b) 


The last fundamental equation which is required for the analysis is 
that which expresses the law of equilibrium of forces acting on a dif- 
ferential wall element in the axial direction. These forces, as shown 
in figure 2(h), lead to the following well-known equation: 


5a , 5 t 


0 


(8) 


Formation of Second-Order Equations 


Equations (7a), (7b), and (8) must he solved simultaneously. 

These three equations contain the four unknowns a, t, I, and w. By 
combining these equations in various ways, three different second-order 
differential equations may he obtained. These three equations eapress 
relations between the displacement I and one of the three vari- 
ables a, T, or w in the following forms: 


5^w E_|_ ^ 

5s^ G as dz " 


(9a) 


a^a 

5s2 


+ 


11 ^ 


+ E’ 


a3| 

5s 5z^ 


= 0 


(9b) 


2 2 
a T E' a T 

as2 G az2 



(9c) 


The first of the above equations may be obtained immediately by 
substituting equations (7a) and (7b) into equation (8) . Equation (9b) 
is most easily derived by differentiating equation (9a) with respect 
to z and multiplying each term by E'. Substitution of equation (7a) 
then gives equation (9b) . In order to derive equation (9c) one may 
differentiate equation (7a) with respect to s and equation (7b) with 
respect to z to obtain 
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or 


and 


_ -pt 


2 

^ ^cr _ ^ ^ y 

E' ^ be 


St _ ^ S V 
^ ~ Ss Sz 


+ G 



Subtracting from the above equation the preceding one gives 

St G So S^I 

Sz E' ^ ~ >,2 

oz 


Differentiating the above equation with respect to z and equation (8) 
with respect to s gives 

S^T _G_ S^g S^g 

Sz2 ■ E* Sz Ss ^ ^^3 

and 

^2 ^2 
G S T G 00 ^ 

E' >. 2 E‘ Ss Sz ^ ^ 

OS 

Adding the above equations and multiplying by E ’ /g gives 
equation (9c). 


The fact that equation (9b) may he obtained from equation (9a^ 
indicates that there is a close mathematical association between the 
axial displacement w and the normal stress a, A stress function F, 
introduced by Hadji-Argyris and Dunne (reference lO), bears a similar 
close relationship to the shearing stress. In extending the analysis 
to semlmonocoque, or stiffened, sections it is advantageous to use 
shear flow as an unknown rather than shearing stress. The shear flow 
and normal-stress flow are defined by the equations 


p = t^ (lOa) 

q = t^T (10b) 

where t^ is the wall thickness and will be considered to be a constant 
for convenience. The stress flows are related to the stress function 
by the equations 
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(11a) 


q = 


Sf 


(ll'b) 


In equation (ll'b) tlie term is the shear flow that may he conqnited 

from the section torque T with the formula from the torsion theory 
of St. Venant 


T 

^ 2A. 


( 12 ) 


where A is the area of the cell formed by the center line of the wall. 
If equations (ll) are substituted into equation (8), the equation of 
equilibrium, it will he seen that the equation is satisfied hy any 
stress function. 


It is now possible to obtain a second-order equation relating F 
to the displacement |. Multiplying equations (?) hy t^ and 
substituting equations (ll) into them gives 


8f 


= E't. 


w 'Sz 


^ ^ 




+ H 


or 


E 

G 


8f 


^ ^ = -E't 


8w 




G 


<3o 


Differentiating the first of these equations with respect to s, the 
last with respect to z, and adding give 


sS' E' aS a^t 


hq 

E* o 
G dz 


( 13 ) 


A second-order equation may he obtained for the shear flow hy differ- 
entiating with respect to z- 


8s‘ 


/8f 

\ E' 1 

^8f _ \ 

■ S 


” ^cy 


8z^ 


= 0 
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The quantity q^ nay he inaerted into the first term since it is a 
function of z only. Substituting equation (lib) gives 







( 14 ) 


It should be noted that this equation is also obt^iined immediately by 
multiplying throu^ equation (9c) by t^. 

It is now seen that the solution to tbe stress problem may be 
obtained by solving any one of tbe four equations (9a), (9h), (13) ^ 
or ( l4) . In each case the first step must be the elimination of the 
tangential displacement ^ . This elimination from equations- (9a) 
or (9h) is obtained throu^^ considerations of equilibrium between 
internal stresses and external loads- On the other hand, the elimina- 
tion of 5 from equations (l3) and (l4) is obtained through considera- 
tion of the necessity for continuity of the axial displacements. Con- 
sequently the process of solution for equations (9a) and (9b) is 
entirely different from that required for equations ( 13 ) and (l4). Thus 
the two possible methods of solution differ, essentially, in the order 
in which one must introduce the concepts of equllibri'um and continuity 
into the analysis. The analysis of Hadji-Argyris and Dunne (refer- 
ence 10 ) involved the solTxtion of equation ( 13 ) while the analysis of 
Von Kannan and Chlen (reference 9), and also that of Duberg (refer- 
ence 18 ), is based i:q)on equation (9a). The analysis of this paper will 
be based on both equations (9a) and (9b). No further consideration 
will be given to equations ( 13 ) and (l4). 


Principal Shear Axes and Associated Section Properties 

In the development of an exact analysis in which axial displace- 
ments, or normal stresses, are considered to be the unknowns, it is 
found to be very convenient to use coordinate axes in the cross section 
which have been called the principal shear axes. Convenient formulas 
for computing the location of these asres were given by Duberg (refer- 
ence 18 ) . At any given point of the wall a tangent to the center line 
may be drawn as shown in figure 3(a). The radius from the origin of 
the principal shear axes to the tangent is indicated as r while the 
angle made by the tangent with the positive direction of the x-axis, 
measured positively in the counterclockwise direction from the axis, 
is Indicated as a. The location of the principal shear axes is defined 
in terms of r and a by the condition that the following three 
integrals must vanish: 


t^ cos a ds = 0 


(15a) 
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sin a ds = 0 


(15b) 



cos a sin a ds 


= 0 


(15c) 


The integration is to he carried over the area of the cross section. In 
the case of a single-cell section such an integration heconies a line 
integral and hence the symbol for a line integral has been ueed for 
convenience. (This integrating symbol would not he aj^ropriate for a 
multicell section.) 

The principal valiaes of three section properties, which are 
associated with the principal shear axes, may he defined hy the 
following formulas: 



(l6a) 

(16b) 

(16c) 


The first of < these section properties was introduced hy Kirste (refer- 
ence 5 ) and called the central moment of Inertia. The magnitude of this 
constant lies intermediate between that of the torsion constant J and 
the polar moment of inertia Ip 


Ip g- le S J 

The constant Ah may he considered as the area resistant to horizontal 
shear and Ay may he considered as the area resistant to vertical 
shear. The stmi of the "horizontal" and "vertical" shear areas is seen, 
from equations (l6h) and (l6c), to he equal to the total area of the 
wall of the heam 


Ag + Ay = A^ 


Aw 



ds 


wi 


where 
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The areas Ag and Ay were also introduced hy Kirste (reference' 5 ) and 
called reduced areas. 


Elimination of Tangential Displacement 

In order to eliminate the tangential displacement I from equa- 
tions ( 9 a) and ( 9 t), it is necessaiy to make further use of the asstmq)- 
tion of rigid hulkheads. Because of this assumption a cross section of 
the beam eijqjeriences "rigid body" displacement under the action of load. 
The horizontal displacement u and the vertical displacement v are 
parallel to the principal shear axes and must in general he regarded as 
the translations of the oi*igin of the principal shear axes. The rota- 
tion, or angle of twist, of the cross section- is indicated hy cp. These 
displacements, and the corresponding contribution that each makes to the 
tangential displacement 5 , are illustrated in figure 3 * Beca\ise of 
rotation about the origin the tangential displacement is given by the 
product rq). Because of horizontal translation of the cross section 
the tangential .displacement is u cos a. Because of vertical trans- 
lation of the cross section the tangential displacement is v sin a. 
These individual effects may be superlnposed to obtain the formula 


S = rep + u cos a + V sin a 


(17) 


If the above equation were used to eliminate 6 , one would obtain 
a relation between the unknowns w or a and the displacements cp, u, 
and V which are also unknown. It is necessary to find a relation 
between cp, u, and v and the applied loads in order to obtain a 
direct relation between w or a and the applied loads. The required 
relations are obtained by equating the section torque and section 
shears, as computed by statics from the external loads, to the section 
torque and section shears as computed by an appropriate integration of 
the internal shearing stresses. (Secondary stresses in indeterminate 
beams will not be discussed.) Thus the section torque T, the hori- 
zontal shear H, and the vertical shear V are related to the internal 
shearing stresses by the follcrwing equations: 



(18a) 



t T COS a ds = H 

W 


(l8b) 



sin a ds = V 


(18c) 
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In equations (18) the shearing stress may he eliminated to obtain 
a relation between the tangential displacement and the loads by 
substituting equation (Tb) . 


G 






ds + G 




ds = T 


G 



cos a 


8w 


ds 



cos a 




ds 


= H 


G 



sin a 


8w 


ds 



sin a 




ds 


= V 


It is now necessary to eliminate the tangential displacement from the 
above equations by substituting equation (l?)* In making this substi- 
tution it will be found that in each equation two integrals vanish 
because of the applicability of equations ( 15 ) which define the loca- 
tion of the principal shear axes. After substitution of equations (l 7 ) 
and (15), the above equations become 

° If 3s + 0 S 3s = T 


G^ty. cos a ^ ds + G cos^a ds = H 


G^ty sin a ^ ds + G 4 ? 


dz 




sin a ds = V 


In the above equations the second integral in each equation is 
seen to be one of the section properties which were previously intro- 
duced by equations ( I6) . Because of the use of principal shear axes 
one may solve the above equations for the derivatives of cp, u, and v 
independently rather than having to solve a simultaneous system. This 
is the reason for employing the principal shear axes as reference axes. 
Inserting the ' section constants and solving for the cross-sectional 
displacements gives 


dcp _ T 
dz “ GIc 



ds 


(19a) 
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IT 


^ = J5_ 

dz ^^2 


— cos a ^ ds 
% J os 


(19b) 


^ e JL 

dz GAy 



sin a 




ds 


(19c) 


The above equations may be substituted into equation (l7) if it is 
first differentiated to obtain 





du 

cos a — + 
dz 


sin a 


dv 

dz 


Siibstituting equations (19) into the above equation gives 


Tr H cos g V sin a 
Sz ^Iq GAy 



cos a 
% 



cos a 


h\r 


ds 


sin a 




sin a 


hw 


ds 


The above equation gives I as a function of the applied loads 
and V, Thus it is seen that this equation may be used to eliminate I 
from equation (9a) bo obtain an equation relating the axial displace- 
ments to the applied loads. Substituting the above equation ?.nto 
equation (9a) gives 


EJ_ 

Ss^ ^ az2 


1 dr 

I ds 
c 





ds 




d cos a 
ds 



cos 



ds - 


1 d sin q jT , . ^ _ T dr H d cos a V d sin a 

Ay ds jP sin a ^ as - - ds ” GA^ ds 


(20k) 

This partial integrodifferential equation, with the appropriate boundary 
conditions, defines the relation between the axial displacements and 
the loads. A corresponding equation for a may be obtained by differ- 
entiating with respect to z and multiplying through by E'. 
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E* 

8s2 G 


1 dr 
Ic ds 



8a 


ds 


% 




cos a 


da 


ds - 


2^ 

Av 


d sin g 
ds 



sin a 


da 


ds 


E* /t dr , ^ d cos a ^ d sin a \ 
G [ic ds Ag 3 s~" Ay ds j 


(201d) 


On tlie rig^t-liand side of equation (20b) t is the torsional load per 
xmit of length, pg is the horizontal conrponent of resultant load per 
unit of length, ana py is the vertical component of resiiltant load 
per unit of length. 


At the root of a cantilever beam, or the end of a fixed- ended 
beam, the boundary condition is that w is zero. At the tip of a 
cantilever beam, or the end of a simply supported beam, the boundary 
condition is that a is zero. This latter condition may be interpreted 
as mpani-ng that dw/dz is zero. A third type of boundary condition 
may arise when shearing or normal stresses of known -distribution ard 
being applied to the free ends of the beam. 


Deflections Due to Shearing Strains 


Before proceeding to a solution of the equation it is worth while 
to draw a few stnq)le conclusions concerning the deflections of the beam 
from the equations of the pre-vious section. Equation (l?) for I may 
be substituted into equation (Tb) for t to ob-bain 

T = G^ + Gr-^ + 6 COB a + G sin a (2l) 

os dz dz dz 


It is con-venient to consider, for the moment, that the beam is loaded 
in such a way as to produce ro-tations only -without translations. 
Equation (21) -would then become 


r 


^ _ T_ dw 
dz ~ G ^ 


dw 


= ^ ^ 


This equation indicates that the angle of -twist may be considered to 
consist of -fcwo parts, the first being d\ie to shearing strains and the 
second due to axial displacements. If this equation is multiplied 
throu^ by -t^r ds and integrated over the cross section, equa- 
tion (I9a) -will be obtained. This indicates immediately that the first 
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term on the right side of equation (l9a) gives the rate of twist 
( change of rotation per unit of length) due to shearing strain, while 
the second term gives the rate of twist due to axial displacements. A 
similar Interpretation can he given to the terms on the right-hand side 
of equations (l9b) and (l9c). 

It is of considerable practical importance to he ahle to compute 
deflections of thin-walled beams due to shearing strains. The proper 
formulas for this purpose are obtained from equations (19) "by setting w 
equal to zero. 


dcp _ 
dz 

T 

GIc 

(22a) 

du _ 
dz 

H 

GAh 

(22b) 

dv _ 
dz “ 

V 

GAy 

(22c) 

These formulas were given by Drymael (reference 8). Equations (22) are 
exact formulas for the .derivatives of the deflections due to shearing 
strain of a beam with rigid diaphragms. It is of considerable interest 
to note that the deflections due to shearing strains are independent of 


the distribution of the sheairLng stress on the cross sections. 

In the case of torsion there is no practical need for computing 
the angle of twist due to shearing strains at various points along the 
span. However, equation (22a) is useful in giving the rate of twist at 
a fixed end where w = 0. This is useful in certain methods of deflec- 
tion calculation. In the case of horizontal or vertical bending, equa- 
tions (22h) and (22c) enable the designer to compute the deflections 
due to shearing strains at various points along the span. This gives 
an indication of the physical significance, as well as the practical 
utility, of the section constants I^,, Ag, and Ay. 

When the axial displacement w is known it may be substituted into 
the Integrals which appear in equations (19) in order to determine the 
deflections. By substituting an approximate solution for w into 
equations (19) one may obtain an approximate solution for the deflec- 
tions. In order to show the relationship between equations (l9b) 
and (19c) and the engineer’s customary equations relating bending 
deflections and loads, the axial displacements may be assimied to have 
planar distribution. It is more convenient to substitute for normal 



20 


KACA TN 2529 


stresses than axial displacements. Equations (l9t) and (l9c) may be 
differentiated with respect to z to obtain 




% 1 r, as ^ 

■ E^J tv cos a ^ ds 

^ t r+ A ^ 

^ sin a ^ ds 


(23a) 


(23h) 


The normal stress in equations (23) may be assumed to be given by 
the flexure formula in the form 


Av ly Ix 


(24) 


In equation (24) the distances x and y are measured from the cen- 
troidal axes of the cross section. The moments Mx and and the 

moments of inertia Ix and ly are coniputed with respect to the cen- 
troidal axes. In general, the centroidal axes, or principal bending 
axes, do not coincide with the principal shear axes. Equation (24) may 
be substituted into equations (23) to obtain approximate formulas for 
the curvatures of the beam. This substitution is explained in detail 
in appendix A and leads to the following equations: 


d^ "Ph ^ ® ^ ® 

^ = E'ly E«I^ 


(25a) 


^2-v- -p^ My. sin 9 cos 9 

^ ^ E-Iy - “ e-I^ 


(25b) 


In these equations 9 is the angiilar difference between the position 
of the principal shear axes and the principal bending axes. In the 
great majority of sections 9 will be small and will have the value of 
zero when the cross section has an axis of synmetry. If 9 is small 
enough that the terms containing sin 9 may be dropped and the cos 9 
replaced by tinity, equations (25) take the following well-known form: 



GAg - E'ly 


(26a) 
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d^v _ ^ 

dz2 “ ■ GAv " E*I^ 


(261) 


In equations (26) the first term on the right-hand side is an exact 
formula for deflections due to shearing strains while the second term 
is an approximate formula for the deflections due to axial displace- 
ments or, as is more commonly stated, due to normal strains. 


Simplifications Due to Symmetry of Cross Section 

A maximum amotmt of simplification arises in the analysis when the 
cross section of the heam is symmetrical about both of the principal 
shear axes. In this case the centroidal axes are also the principal 
shear axes. The general integrodifferential equation for w may be 
replaced by three independent equations in which the actions of torsion, 
horizontal bending, and vertical bending become completely separated. 
From considerations of double symmetry it is found that equation (20a) 
must be replaced by the following three equations: 


^^w E* ^^w 1 dr ^ T dr / „ v 

;^„2 G ^^2 ' Ic ^ ~ ' GIc de 


S^w . E' d^w 1 d cos a 




® az2 Ag ds 




t^ cos a 


Bw 


ds = 


H d cos g 
GAg ds 


(27b) 


5^w ^ ^ 5^w 1 d sin g 

Ss2 G ^^2 Ay ds 




V d sin a 
GAy 


ds 


(27c) 


Equation (27a) was derived and solved by Von K^rrn^n and Chien (refer- 
ence 9)- The development of equations (27) is given in appendix B. A 
corresponding system of equations may be derived for the stress distri- 
bution in a beam with doubly symmetrical cross sections. These equa- 
tions are as follows: 


^2(7 ^ S2(j 

3 s 2 G ^2.2 


1 dr T, j E't dr 


GI^ ds 


(28a) 
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0 3^2 


_l_ \_WO VJU Wfv , 

5;— 


COS a 


^cr 


ds 


E’p 

H d cos g 

GA,g ds 


(28b) 


p p 

8 a E ’ 8 a 

8s^ ^ 8z^ 




d sin g 
ds 



sin a 


8a 


ds 


^ ^ d sin g 
GAy ds 


(28c) 


Every 'term, of equation (27a) or (28a) is anti symmetrical about both the 
X-axis and the y-axis. Every term of equation (27b) or (28b) is sym- 
metrical about the x-axis and anti symmetrical about the y-axis. Every 
term of equation (27c) or (28c) is symmetrical about the y-axis and 
anti symmetrical about the x-axis. 

If a cross section has only one axis of symmetry, the general 
equation for v is replaced by two equations. For the case of a 
cross section which is symmetrical about the x-axis the equations take 
the following form: 


p p 

8 w E' 8 w 

3s 2 0 a^2 


1 d cos 




ds 




cos 


8w 


ds 


H d cos g 
GAg ds 


(29a) 


8^w E* 8^w 

8s2 ^ 8z2 


1_ 

I 


dr 

ds 




8w 


ds 


Ay 


d sin g 
ds 



sin a 


8w 


ds = 


T dr V d sin a 

GIq ds GAy ds 


.(29b) 


Equations (29) are derived in appendix B. Equation (29a) indicates 
that the solution for axial displacements due to horizontal loads may 
be considered separately from the action of vertical and torsional 
loads. Every term of equation (29a) is symmetrical about the x-axis 
while every term of equation (29b) is antisymmetrlcal about the x-axis. 

For the case of a cross section which is symmetrical only about 
the y-axis, the equations for w take the following form: 


2 2 
8 w ^ 8 w 

8 s^ ^ 8z^ 


1 d sin g jC 
Ay ds J 


Bloa|2 dB 


V d sin a 


ds 


ds 


(30a) 
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Ss2 G Sz2 






ds 




d 


cos g 
ds 



ty cos a 


hv 


ds = 


T dr _ H d cos a (SOt) 

GIj. ds ~ GAg ds 

Eguations (30) are derived in appendix B. Every term of equation (30a) 
is symmetrical about tbe y-axis while every term of equation (30h) is 
anti symmetrical about the y-axis. It is obvious that equations for the 
normal stress, corresponding to equations (29) and (30), maybe readily 
derived. 

A few comments concerning the effects of symmetry upon deflections 
may be of some interest. Considering a cross section having double 
symmetry, the solution for w from equations (27) may be substituted 
into equations (19) to obtain the deflections of the beam. It will be 
found that torsional loads produce rotations only, horizontal loads 
produce horizontal translations only, and vertical loads produce verti- 
cal translations only. For the case of a cross section which is sym- 
metrical about the x-axls, the solution of equations (29) may be sub- 
stituted into equations (19) as before. It will be fotind that hori- 
zontal loads produce only horizontal displacements. However, vertical 
loads produce both vertical displacements and rotations, and torsional 
loads also produce both vertical displacements and rotations. In the 
case of a cross section which is symmetrical about the y-axis, vertical 
loads produce only vertical displacements but there is an interaction 
effect between horizontal loads and torsional loads. Such interaction 
effects do not occur in elementary beam theory wherein one considers 
only primary stresses. 


Solution of Homogeneous Equation 

The solution of the nonhomogeneous equation which governs the 
axial displacements, or the normal stresses, may be determined by 
assuming that the unknown may be e3q)anded into an infinite series of 
appropriate orthogonal functions which individually satisfy the 
boTindary conditions and the conditions of continuity. The orthogonal 
functions which are appropriate for the analysis imy be determined by 
considering the solution of the homogeneous equation and the boundary 
conditions. In the case of cantilever beams, simple beams, or fixed- 
ended beams, the end conditions are homogeneous. The transverse dis- 
tribution of the solution must be continuous and single-valued, or 
periodic. 

The case which is most easily understood and capable of being given 
the simplest explanation of the method of solution is that of a beam 
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with douhly synnnetrlcal cross section acted t^jon hy torsional loading 
only. It was this case which was treated hy Von Karman and Chien, 
although their considerations were limited to a cantilever heam of 
semi -infinite length acted upon hy a torgue at the tip. In the present 
analysis the loading distribution and end conditions will remain 
arbitrary. 


The homogeneous equation governing the axial displacements in a 
cross section with double symmetry is obtained by setting the right- 
hand side of equation (27a) equal to zero. 


E' S^w 1 dTvdCj. ^ 


The corresponding equation for normal stresses becomes 


(31a) 



^2 

E* B a 



1 dr 
Ic ds 



Bn 


ds 


0 


(31t) 


A standard well-known method for the solution of homogeneous partial 
differential equations is to assume that the unknown may be e^ressed 
by a separation of the variables. The same method may be employed in 
solving the above integrodifferential equations. The displacement and 
stress may be assumed to have the following foims: 


w = f(s)g(z) (32a) 

a = f(s)h(z) (32b) 


where 


h = E' 


is 

dz 


(32c) 


The relationship between g and h arises from the relationship 
between w and a as ejqpressed by equation (Ta) • It is also clear. 
from equation (7a) that, if w and a can be e3q>ressed by a separa- 
tion of the variables, then the transverse distribution of both func- 
tions must be the same. 

Substituting equation (32a) into equation (31a) gives 


g 


. o ^ a -Q T_Sar 
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Dividing throng by the product fg and transposing the first and 
third terms gives 


11 




di 


(33) 


Since the left side is a function of z only while the right side is 
a function of s only, each side may be equated to a constant inde- 
pendently. Equating the left side to gives 


or 


where 


2 2 
^-^g = 0 

,2 E' ® 
dz 


- ."g = 0 

dz'^ 


J 





(3^Ti) 


The corresponding equation for the spanwise distribution of a is 
obviously 

- A • 0 (31k) 

Equating the right-hand side of equation (33) "to gives the 

following equation for the transverse distribution of w or cr: 



X^f 


1 dr J? df 


(35) 


It is convenient for the analysis to introduce the following integration 
by parts; 




tvff ds 


(36) 


The bracketed term, which usually arises in an integration by parts, 
vanishes here because f must be continuous and periodic. Althoti^ 
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eqiiation ( 35 ) is an integrodlfferential equation it is convenient in 
certain phases of the analysis to consider the integral in the eqmtion 
as a section constant for the heam. If the stress distribution were 
known it cotild he substituted into the integral which, in turn, could 
be evaluated to obtain a scalar constant that would be dependent only 
upon the properties of the cross section. This section constant, which 
is associated with torsional loading, may be indicated as Ct and 
defined as follows: 


Using equation ( 37 ) 
following forms: 



, equation ( 35 ) converted into either of the 


d^f 

ds2 


X^f 


Ic ds J 




(38a) 



'"T dr 

X f + ■= — — — = 0 
I ds 
c 


(38b) 


It is of some interest and benefit to draw a comparison, at this 
point, between the present problem of detennining transverse stress, or 
displacement, distributions, and the natural vibrations of a thin ring. 
The analogy is of a physical nat-ure only, since the problems do not 
correspond mathematically. The ring that is to be considered lies in a 
plane and has a doubly symmetrical shape. It is to be considered as a 
free body vibrating in space without a gravitational field. The only 
vibrational displacements that are to be considered are those that are 
normal to the plane of the ring. There will be four infinite sets of 
modes that could be excited as follows: (a) Doubly symmetrical, 

(b) doubly anti symmetrical, (c) symmetrical about the x-axls only, and 
(d) symmetrical about the y-axls only, where the x-axis and y-axls are 
in the plane of the ring. In a similar manner there will arise four 
infinite sets of nattiral "displacement modes" or "stress modes" in a 
cross section having double symmetry. Each type of mode will be 
"excited" by a particular type of loading. The relationship between 
type of mode and loading is, as follows (see appendix B) : 
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Mode 

Loading 

(a) 

Doubly symmetrical 

y Axial force 

(b) 

Doubly anti symmetrical 

y Torsion 

(c) 

Symmetrical about x-axis — 

y Horizontal bending 

(a) 

Symmetrical about y-axis — 

> Vertical bending 


The analogy with a vibrating ring is useful in suggesting terminology, 
in visualizing the stress distributions that are being determined, and 
in suggesting calculation techniques for idealized stiffened cross 
sections that are to be treated in the second section of the paper. 

Solutions of equation ( 38 a) exist only for certain characteristic 
values of the parameter X. Ordinarily characteristic nunibers are 
determined by boundary equations. In the present case the function f(s) 
is governed by contintiity and periodicity conditions which do not 
provide a characteristic equation, (in the case of a section with 
comers, the condition of continuity at the comers provides the charac- 
teristic equation.) The characteristic equation was obtained in a 
simple manner by Von Karman and Chien by setting the line integral in 
equation ( 38 a) equal to unity. The resulting equation determines an 
infinite set of values of X. Inspection of equation ( 38 a) shows that, 
if the line integral is replaced by any convenient constant, then the 
eqiiation becomes a nonhomogeneous differential equation which can be 
readily solved. The solution for f(s) will contain the number X 
and hence there will be an infinite number of stress modes. 

UnfortTonately the simple method used by Von Karman and Chien to 
form the characteristic equation cannot be extended directly to arbi- 
trary cross sections. Consequently a method of forming the character- 
istic equation will be presented which can be extended to arbitrary 
cross sections. Equation ( 38 b) suggests the following form of the 
general solution; 


f = A sin Xs + B cos Xs + CijCp(s,X) 


(39) 


where cp remains to be determined. Substituting equation (39) into 
equation ( 38 b) gives the following equation for qp: 




1 dr 

Ic 


(iK)) 


The function cp is the particular integral of equation (4o) and is a 
function of the coordinate s and the parameter X. 
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It is now necessary to show that the constants of integration A 
and- B must vanish for noncircular sections. The stress mode f is 
douhly anti symmetrical. It is convenient to choose the origin of the 
coordinate s on one of the axes of symmetry. The term sin X.s may 
he made douhly anti symmetrical hy a proper choice of values of X. At 
first thought this would seem to he an apprcroriate method of determining 
the characteristic values of X. Since dr/ds is anti symmetrical it 
is seen from eqmtion (to) that <p must he anti symmetrical. However, 
the term cos Xs is symmetrical about the axis upon which the origin 
of s is chosen and, hence, B must he set equal to zero. 

In order that the term sin Xs shall he anti symmetrical it is 
necessary that X have one of the values X^ given hy the following 
formula: 


X^ = ^ n = 1, 2, 3, . . • (^1) 

where L is the circumference of the center line of the wall of the, 
cross section. The geometric quantity dr/ds may he assumed to he 
expanded into a Fourier series as follows: 


dr \ — . torts 


(te) 


The function cp(X, s) may also he expressed as a Fourier series 


'Pn(V^) =^^iW ^^3) 

Soibstituting equations (to) and (to) into equation (to) leads to the 
following formula for hi: 


^1 



where 


m 



i = 1, 2, 3, 


(to) 


For each value of cp^, as defined hy equation (to)^ there will he one 
term in the series for which ooi equals X^ and hence for which hi 
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becomes infinite. Hence tbe only possibility for a finite solution 
for qpn to exist is for all of the aj[ values to be zero. TMs means 
that dr/ds must vanish, which would be true only for a circular sec- 
tion. An examination of equation (to) shows that if dr/ds vanishes 
the function cp must either be zero or be proportional to the quan- 
tity sin Xs. In the latter case q) would add nothing to the formula 
for f as given by equation ( 39 ) and hence may be assumed to be zero 
for a circular cross section. For the noncircular section it is iB5>os- 
sible to obtain a finite solution for cp by chCKJsing X to maloe 
sin Xs antisymmetrical. The only alternative in this case is to choose 
A' equal to zero. Thus the following formulas for the stress distri- 
bution have been derived: 

Circular section: 


f = A sin Xs (toa) 

Noncirciilar section: 

f = C^^(X,s) (tob) 

In the case of the circular section X is given by equation ( 4 l). For 
the noncircular section a characteristic equation must now be derived. 

Multiplying through equation (tob) by dr/ds and integrating over 
the cross section give 




ds 


or 




or 

VP £ j Ct = 0 (to) 

Equation (to) te regarded as a homogeneous algebraic equation in 
which Ct is the unknown and thus remains arbitrary. The term Ct 
plays the role of an anqolitude factor which must be determined subse- 
quently by introducing a normalizing condition for the stress modes. 

The coefficient of Ct must vanish and hence 
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^VP If ^ 

Equation (it8) is tlie characteristic equation which determines the -values 
of X and agrees -v/ith that given hy Von Karman and Chien. The aho-ve 
deri-vation of the characteristic equation can he extended to arbitrary 
cross sections as will he shown in the next section. 


The section cons-tant is determined hy introducing the following 
normlizing condition: 



( 1 ^ 9 ) 


!This normalizing condition is more conTCnient in the analysis, for hoth 
monocoque and semimonocoque sections, than the engineering types of 
normalizing laws -which might he s\iggested hy the -vibration analogy. The 
method of solution for the na-tural stress modes should now he clear. 

The function q) is detennined in terms of X hy sol-ving equation (4o) . 
Eacpansions in Fourier series should he -useful for this purpose. The 
functions cp and dr/ds are then suhsti-tu-ted into equation (W) to 
determine the characteristic -values of X. The stress modes are then 
known except for the coefficients -which are determined to satisfy 

equation ( 49) • 


Stress Modes for Arbitrary Cross Sections 

The solution for the stress modes and characteristic numbers for 
arbitrary cross sections will he stated briefly. If the ri^t-hand 
side of equations (20) is replaced hy zero and a separation of -vari- 
ables is introduced, the folio-wing equation -will he found to govern the 
trans-verse distribution of displacements or stresses: 

£1 ^ ^ ^ d cos g S d sin g _ 

ds2 Ic % ds Ay ds 



-where 


Ct 



(50h) 


Ch = 



d cos 


ds 


a 


ds 


(50c) 
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Cy =;f ^ — "d - 3 ds (50d) 

The solution of equation (50a) may he written as follows: 

f = CrjTcpj(l,s) + Cg<j^(l,s) + Cy<Ty(1,s) (51a) 

where (Pg, and cpy^ aire governed hy the following equations: 


. -\2_ Id cos a 

" 'T' A 


ds" 


%--i2 


ds 


(51c) 



2 

X cpy 


1 d sin g 
Ay ds 


(51d) 


Terms containing constants of integration in the formula for f have 
been omitted since they may be shown to vanish for noncircular sections 
by arguments similar to those used in the previous section. Since, in 
the present case, bending as well as torsional action is being con- 
sidered, the complete solution must include the planar stress distri- 
butions of ordinary engineering beam theory. These primary modes must 
satisfy equation ( 50a) . The plaimr distribution may include the effect 
of an axial force as well as, bending moments. 


If a planar distidbution of stress is assumed, it will be found 
that the section constant Ciji has the value zero by evalxiating the 
integral which it represents. The value of X which is associated 
with the planar modes is also zero. If cpw and cpy are defined to 
have the linear distributions a + x and b + y, respectively, in 
equation ( 51 a) for f, this planar distribution will be found to satisfy 
equation (50c) if X and are set equal to zero. This existence 

of a solution for a zero value of X is a featirre of the present 
problem which is not found in the solution of homogeneous second-order 
differential equations with homogeneous boundary conditions (the 
Sturm-Liouville problem) . 


The characteristic equation iisy now be developed by the same 
method that was used in the previous section. Equation (51a) must be 
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multiplied tliro\igh successively "by the following quantities: 
dr/ds, d cos a/ds, and d sin a/ds. Integration over the cross sec- 
tion in each case gives the following three linear algebraic equations: 


Ct = t^ g da + VPh If ds + VPv ^ (52a) 






d cos CL _ d cos a 

d^ + Ch 


^ Vte - - g - ■!= + W 


d cos a 


ds 


ds 


( 52 b) 


^ ^ .X j. <3. sin a , ^ .X j. ^ sin a , . r, P ^ <3 sin ct , 

Cy = Crpy tj/Pr ds + Cgy t^ ds + Cy^ t^/py ^ ds 

( 52 c) 


Equations (52). are a homogeneous system of three algebraic equations in 
which Cg, and Cy are the unknowns. The integrals which appear 

are functions of X and may be regarded as scalar coefficients of the 
equations. After transposing all terms to one side of the equations 
the determinant of the coefficients may be isolated and set equal to 
zero. This provides the characteristic equation which determines the 
values of X. The calculation task of determining the characteristic 
values of X is clearly a formidable one for arbitrary cross sections. 
For each value of X a solution of equations (52) exists for relative 
values of Cj, Ch, and Cy. Their specific magnitudes must be deter- 
mined from the normlizing condition as ejcpressed by equation (^ 9 )* 


Orthogonality of Stress Modes 

In order to obtain a sol\rtion to the general nonhomogeneous equa- 
tion for stresses or displacements it is necessary to determine the 
orthogonality properties of the stress modes which have been defined. 
Since this argument follows customary well-known mathematical methods, 
it will be given only for the doubly anti symmetrical modes (the torsion 
case) in a beam having a cross section with double symmetry. 

Equation ( 38 b) may be written for the mth and nth stress modes as 
follows : 


d^f 


ds‘ 


r + Cfj 


m 


+. 


^Tm dr 


(53a) 
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n 


ds^ 


^Tn dr 


S = 0 


(53b) 


Equation (53a) is to be multiplied throu^ by equation (53b) by fjjj, 

and botb equations are to be integrated over the cross section. 

P d^f_ p P 

= 0 (5l*a) 


P 2 ^Tn^Thn 

y fm ^ ^ V^t^ ds + = 0 (54b) 


The first term in each of these equations may be integrated by parts 
to obtain 




df. 


m 


ds 


t^ ds 



ds 






The bracketed tenns vanish since the stress mode and its derivative 
must be continuous and periodic. These formulas show that the first 
terms of equations (54a) and (54b) are identically equal. Hence, by 
subtracting equation (54b) from (54a), one obtains 


■3= - 0 

Since the X’s form' a discrete set of ntmibers, the first factor in 
the above equation does not -vanish and the following orthogonality 
condition is ob-fcained; 

^Vn-tw^^=0 (55) 

A similar de-velopment of equation (55) can be formulated for arbitrary 
cross sections. 

In deri-ving the solution of the nonhomogeneoua equation for 
stresses or displacements it is convenient to ha-ve one additional 
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relationsMp wMcli may now be derived. Ifelcing use of the normalizing 
condition as given hy equation (l<-9) may be seen that, when m = n, 

the second term of equation (5^b) is equal to or . When m 

does not equal n, the second term of equation (5^b) vanishes because 
of equation (55) • Thus equation (5^b) may be rewritten in the following 
manner: 



2 

d f 


n 


ds^ 


tw ds 


^Tm^Tn 

Ic 


"\i ^mn 


(56a) 


where 



(56b) 


Solution of Wonhomogeneous Equation 

Again the action of torsional loads on a beam having cross sec- 
tions with double symmetry will be considered first. The equations 
that must be solved are equations (27a) and (28a) . The displacements 
and stresses may be assumed in series form as follows; 


w = ^ fn(s)g^(z) 
n=l 


(5Ta) 


00 

CT = 51 


n=l 


Substituting equation (5?a) into eqmtion (27a) gives 




L 


Ic ds 


(5Tb) 




T dr 


w 


ds 


Interchanging the order of summation and integration in the third term 
of this equation and introducing as defined by equation (37) give 






n 


T dr 


GIc ds 
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The preceding eguatlon must now he multiplied, throu^ hy and. 
integrated, over the cross section. 


z-f 


. th 

^ 2 
de^ 


t^ ds + 




i ^ ClnSn - - lin S tw 


The integrals which occur in the second series of this equation may he 
evaluated from the normalizing condition and the orthogonality condi- 
tion for the stress modes. All terms of the series vanish except the 
mth term. After introducing Cij^ and rearranging terms the equation 

becomes 


11 ^ ^ 
® dz^ 




"b^ ds + 


^Tm^Tnj 




The qmntlty in the parentheses has a finite value only when m = n. 
Substituting equation (56) into the above equation will give the equa- 
tion governing the spanwise variation -of the mth displacement mode. 

The subscript m may be changed to n for convenience to give the 
equation for the nth mode as follows: 

E' d^gn , 2_ CtuT 

G dz2 ” ^ = ■ GIc 


Multiplying through by G/E’ and introducing ^ as defined by equa- 
tion ( 3 ^b) gives 


^ 8n 2 _ Gj^T 

^^2 ' - E'l^ 


(58a) 


The coire spending equation for the stress distribution is 



Mn ^n = 


C^nt 


(58b) 


The complete solution can now be obtained. After computing the 
transverse stress modes as previously explained the constants Cijq 
and Pn are known for each mode. Equation (58a) may tfien be solved 
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for the spanwise variation of the displacements due to the applied loads. 
The solution of equations ( 58 a) or ( 58 b) may he -written in either of 
the following forms: 


= Ae^^ + Be + tKs) (59a) 

= A sinh |j.^s + B cosh p^s + ilfCs) (59t) 

The coefficients A and B nay he determined from the boundary condi- 
tions after the particular integral -<!f(s) has been determined for the 
gi-ven loading. 

In the case of a heam of arbitrary cross section without symmetry 
the same mathematical method may he applied to equations (20) . The 
resulting equation for the spanwise -variation of axial displacements is 

^ 2 _ ^Tn^ _ 

^^2 “ n ’ E-A^ " E'A^ 

A corresponding equation may be -written for the stress distrib-ution. 
Equations (59) also gi-ve the form of solution for equation ( 60 ) . It is 
of interest to note that a zero--warping axis for the nth mode may be 
located by setting the ri^t-hand side of equation ( 60 ) equal to zero. 
Howe-ver, it is not apparent that there is any practical value to be 
gained by locating this axis in the present method of analysis. 

Associated with each normal stress mode there is also a shearing 
stress distribution. These sheeiring stresses are computed for each 
normal stress mode independently by integrating equation ( 8 )^ the equa- 
tion of equilibrium. This integration is also required in primary 
stress analysis and hence is well-known. The constant of integration 
for the secondary shear flow is determined from the condition that the 
resul-fcant torque must be zero. 



H - BEAMS HAVING TEEN-WALLED STIFFENED CROSS SECTIONS 

Introduction 


In section II a method of analysis is to be presented which is 
directly applicable to beams of practical cross sections such as occur 
in aiip)lane wings and fuselages. The beam -will be assimied to have a 
thin--walled stiffened cross section which may be of single-cell or 
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raulticell design. Tlie type of cross section "wiiicTi Is actually to be 
anal3Tzed is the idealized stiffened section wMcb results from the 
customary engineering idealization that is used in primary stress 
analysis of airplanes. The loading distribution and end conditions 
will again be considered as arbitrary. 

As before the beam will be assumed to contain closely spaced 
internal rigid diaphragms which stabilize the shell and preserve the 
cross-sectional shape. There are two sources of error in the present 
analysis. The first is the idealization of the cross section and the 
second is the assumption of rigid bulkheads. As explained in section 1 , 
it is possible to obtain "approxlamte” solutions or "accurate" solutions 
for stresses and displacements. Only accurate solutions will be 
considered herein. 

The introduction of generalized difference equations permits the 
analysis of a semlmonocoque, or stiffened, section to be developed in a 
manner that is exactly analogous to the method of solution which has 
been previously given in te 3 rms of analytic functions for monocoque 
sections. The use of generalized difference equations provides certain 
conveniences and slrplifications which are not found in the customry 
algebraic treatments of the stiffened beam. The differencing symbols 
have been specifically defined for application to the stress -analysis 
problem of idealized multicell cross sections. 


Symbols 

aj area of jth stiffener 

A^ central area associated with tth web 

Ag total area of stiffeners 

j stiffener number 

k web number 

length of kth web 

Pj axial force in jth stiffener 

warping stiffness of kth web 

geometric discontintilties at jth stiffener 

[a ] matrix of stiffener areas 
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[B] matrix of -warping stiffnesses 

[r] matrix of geometric discontinuities 


Difference and Summation Definitions 

Tlie idealized semimonocoque section -which is to he analyzed has a 
finite number of stiffeners -which carry all of the normal stress and 
the shear flows are distrlhuted uniformly o-ver the connecting wehs. 

The idealization is illustrated in figure 4 where a small portion of 
the cross section of a heam is shown. The flange and stringers of the 
true section are replaced hy idealized stiffeners centered on the -wehs. 
The idealized stiffeners also include an "effective -width" of -weh 
material as is commonly employed in primary stress analysis. In the 
idealized section it is asstnned that the -wehs -which connect the stiff- 
eners are strai^t. The curvature of the -wehs is neglected. The 
Idealization of a cross section for secondary stress analysis need not 
he highly accurate for design purposes although for research studies 
some Improvement in the acctiracy of the idealization may he found 
desirable. 

The normal stresses on a cross section act over the eirea of the 
stiffeners only. The normal stress acting on the jth stiffener, 

may he multiplied hy the stiffener area aj to oh-tain the axial 
force pj on the stiffener. The axial forces acting on three stiff- 
eners on a strai^t wall are shown in figure 5(a)* The shearing 
stresses are distributed -uniformly along each individ\aal weh. The 
shearing stress acting on the kfch -weh, may he multiplied hy the 

-wall thickness t-fj; to oh-bain the shear flow acting on the weh. 

The distribution of shear flow over a cross section is gi-ven hy a 
stepped diagram as ill-ustrated in figure 5(h) . 

The aho-ve definitions in regard to stresses are equally izseful in 
primary or secondary stress analysis. In the secondary stress problem 
displacements must also he considered. Because of the fact that each 
idealized -weh has been assumed to he straight, the -tangential displace- 
ment of points on the fcth weh -will he a cons-tant for the weh and 

the trans-verse distrihirtion of -will he a stepped diagram. The 

function -will thus behave in the analysis in a manner similar to 

the stress and the product may he regarded as a "flow-type" 

quantity. The axial displacements must he considered to he defined at 
the stiffeners, the displacement wj occurring at the jth stiffener. 

It must he assumed that the axial displacements have a linesir -variation 
over each weh. The resulting distrihiitlon of axial displacements is 
gi-ven hy a hroken-line diagram as shown in figure 5(c). The distribu- 
tion of axial displacements must he a continuous function from physical 
considerations . 
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All guantities entering into the analysis which are dependent upon 
the transverse coordinates are completely defined hy giving their 
values either at the stiffeners or at one point on each weh. Those 
functions which are defined over the wehs are uniformly distributed 
over each weh. Those guantities which are defined hy their values at 
the stiffeners are given a subscript Those gmntities which are 

defined hy their values over the wehs sire given a sxibscript k. 

Just as two types of functions arise in the analysis, so it is 
necessary to define two types of differencing symbols. The symbol Aj 
indicates a differencing operation which produces a function that is 
defined hy its values at the stiffeners. The symhol Afe indicates a 
differencing operation which produces a function that is defined over 
the wehs. The differencing symbols may he defined most conveniently 
hy applying them to particular functions. The difference A^wj may he 
defined as the difference between the value of wj at the forward 
stiffener and the value at the rearward stiffener relative to the 
tth weh. If the forward stiffener is stiffener a and the rearward 
stiffener is stiffener h, the difference A]jWj is defined hy 

= Wa - Wb (61) 

This definition is illustrated in figure 5 (c). The above difference 
must he divided hy the length Lt of the kth weh to obtain the slope, 
or derivative, of the function. The correspondence between derivatives 
and differences is as follows: 

Throughout the analysis for multicell sections the positive direction 
for the coordinate s and all "flow-type” gueuitities is counterclock- 
wise along the outer shell and in the positive direction of the 
coordinate axes on the interior wehs. 

The symbol Aj is defined as the difference between outflow and 
Inflow at the jth stiffener. This definition of a differencing symbol 
is of a rather general nature since it may he applied at a flaige 
having any nimiber of connectiig wehs. The difference Ajgjj. is illus- 
trated in figure 5 (h) for a stringer with two connecting wehs. If the 
weh ahead of the Jth stiffener is weh c and the weh behind the Jth 
stiffener is weh d, the difference Ajg^ is given hy 

'’c - % 




(62a) 
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In this case no clear coirrespondence between derivatives and differ- 
ences is apparent. A second example of the difference Ajq-[j. is shown 


in figure 6 where the jth stiffener is connected to webs a, 
The positive directions of the shear flows are indicated, 
ence A.g, becomes 


b, and c. 
The differ- 


= % - '3a - 

In addition it is necessary to introduce summations for the 
idealized stiffened sections corresponding to the variotis integrals 
which occurred in the analysis of section I. Two types of stoumations 
must be introduced to indicate summation over all of the stiffener 
areas or summation over all of the web areas. Three simple examples of 
summations over the stiffener areas are as follows: 




(63a) 


<J J 


(63b) 


- Y1 = II pjyj 

d d 


(63c) 


It is apparent that these exumnations give the axial force and bending 
moments on the cross section. Three simple examples of the second type 
of summation are as follows: 


T 


- Y1 VkVk = II ’kVi 


( 64 a) 


H = m cos = 2Z cos 

k k 


(64b) 


V = ZZ Vkht '^k = I— ‘3kh: 

k k 
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These summations give the section torque and shears. The geometric 
quantities cos and sin a-jj. are constant over the tth weh 

since the web is assumed to be strai^t in the idealized section. 

It may be worth while to note that the introduction of generalized 
difference equations in one-dimensional problems would be of no practi- 
cal value. In such problems difference equations become useful only 
when the regularity of the structural layout leads to "standard differ- 
ence equations" such as those treated in chapter XII of reference 19- 
In problems of two or more dimensions, however, it is found that gener- 
alized difference equations are of considerable value. They not only 
provide a mathematical method of complete generality for deriving the 
equations governing the axial displacements or stresses from the basic 
equations of elasticity but, more importantly, they make clear the 
possibility of using the method of separation of variables in spite of 
complete irregularity of the structtiral layout which may occur. .The 
method of analysis for an Idealized stiffened multicell section may be 
carried out, step by step, in a manner exactly parallel to the method 
used in section 1 for the single-cell monocoque section which is more 
easily understood. The method of separation of variables leads auto- 
matically to the natural stress modes for the cross section. The pos- 
sible utility of generalized difference eqmtions in other planar or 
spatial structures, such as grldworks, remains to be investigated. 


Fundamental Equations of Elasticity 


Since the present analysis is exactly parallel to that given in 
section I, it need be stated only briefly. The analysis will be made 
clearer by freqiient comparison with corresponding equations in sec- 
tion I. The stresses are related to the strains by the formulas 



(65a) 

'^k = ^’'k 

(65b) 


The strains are related to the displacements by the following formulas: 


e 




(66a) 




(66b) 


Equations (66) are comparable with equations (6) of section I. From 



h2 


NACA TW 2529 


equations ( 65 ) and (66) it is seen that the stresses 
in terms of the displacements hy the formulas 

Oj = E 



nay he expressed 

(6?a) 


= — A. w. + G ^ 


k - Lfc Vj 


The stress flows are given hy the following formulas: 




(6Th) 

(68a) 


Gt-jj S 

^k = ■‘^k'^k = + ^■‘^k ^ 


= Gp^wj 


^k 


(68b) 


where 


Pk 


Lk 


(68c) 


The quantity GPfc nay he called tte warping stiffness of the kth weh. 
The forces which act on a differe^ial length of a stiffener are shown 
in figure 7- The equation of equilibrium for the stiffener hecomes 


dp. 

IT ° 


This equation corresponds to equation (8) of section I. 


( 69 ) 


Developnant of Equation for Axial Displacements 

In the present analysis the equation will he derived only for the 
axial displacement. Equations for the other functions may he readily 
derived. Substituting equations (68) into equation ( 69 ) gives the 
following equation relating the displacements. 
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Ea^ ^) = ° 

This equation corresponds to equation (9a)- The second term, in this 
equation is a second-order difference quantity. 

Before eliminating the tangential displacement it is necessary to 
introduce the definition of the principal shear axes and the associated 
section properties. The principal shear axes are located hy the condi- 
tion that the following three sunmations must vanish: 

y~ “k = ® 

k 


5Z\^k^ sin = 0 (Tlh) 

k 


y tkLjj cos sin *s 0 (Tic) 

k 

The associated section properties are as follows: 

t 

Ic = ^ t^jLkrk 

k 


% = ZI COS^Ojj. 
k 


\ - II W B±n\ 

Convenient methods for calculating these section properties are given 
in appendix C. 

Corresponding to equation (17) of section I the tangential displace- 
ment of all points on the kth weh is related to the cross-sectional 
displacements hy the formula 


(72a) 

(72b) 

(72c) 
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+ u cos ojj. + V sin a^j,. (73) 

The section torque and shears have been ej^re seed in terms of the shear 
flows hy equations (64). If equation ( 68 b) for the shear flw is sub- 
stitiited into equations ( 61 ^), the following equations are obtained: 


COB + G / “h ' 5 z“ 

k ^ 


G 21 "^k + G 21 -tk^c sin ai, ^ = V 

k ^ 

Substituting equation (T 3 ) and omitting those summations ^ch vanish 
according to equations (Tl) , the above equations become 

2Z ® S ^ ^ 

k K 




G 2 Z t^ ® Iz XZ '’^k^: ^ 

k t; 

Introducing the section properties as defined by eqmtions (72) gives, 
after rearranging terms, 


dz GIc Ic ^ ^ 


^ = -1 ~ S tv cos Ctv^W^ 

dz GAgr Ah ^ J 
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i^5 


^ = _L. - i 

dz GAy Ay 




(74c) 


These eguations were given hy Drynael (reference 8) . 

The required formula for Ijj. which will permit its elimination 
from equation (70) may now he obtained by differentiating equation (73) 
with respect to z and then substituting equations (74). 

81^ Tri, H cos ai- V sin ai- ri, — 


cos — sin ajj. — 

2 _ -tfe cos a^wj - — ^ }_ t^ sin a^w^ 
n k '' k 

This formula may be siibstituted into equation ( 70 ) to obtain the 
following equation governing w^: 

d^w 

+ I ■ T - "h ^JPHk - 

QZ ^ k n k 


Ay ^j^Tk ^ fVk'Vj “ “GIc " GA2 " GAy ^j^Vk 


In equation (75) the flow-type quantities and have 

been introduced in accord with the following definitions: 


%k “ "^k^k 

( 76 a) 

Phj^ = t^ cos 

( 76 b) 

= t^ sin 

( 76 c) 

Eqmtion (75)^ which is for arbitrary 
replaced by simplified groups of equations 
or double symmetry just as in section I. 

cross sections, may be 
for sections having single 


In section I an integration by parts was introduced after the 
variables had been separated. This integration by parts could have 
been introduced earlier in the original partial integrodifferentlal 
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equation "but there -was no apparent reason for so doing. In the present 
case an appreciable simplification of the form of the equation may he 
obtained by performing a summation by parts (corresponding to an 
integration by parts) providing one simultaneously introduces new 
symbols for such geometric differences as A^p^. Consequently the 

following definitions are introduced: 


^Tj - "^jPTk = ^3 f^k^k) 

(77a) 

yR3 = ^.iRSk = ^d(^^k “k) 

(77b) 

= ^JPvk = “k) 

(77c) 


Convenient methods for computing these geometric discontinuities are 
given in appendix D. The operation of summation by parts may now he 
performed on the summations which occur in equation (75) to obtain the 
following formulas: 


= - XZ'^J^jPTk = ' (T8a) 



The validity of the above srmimations by parts is demonstrated in 
appendix E. Equation (75) nay now be written in the following form: 













T 

Gflj, 


H V 

GAg ’hJ " GAy 


(79) 
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Simplifications Due to Symmetry of Cross Section 


Tbe arguments wMch permit the replacement -of the general equation 
for axial displacements hy groups of siu5)lified equations are identically 
the same as those used in section I. Hence it is not necessary to 
repeat the eirguments hut merely to state the equations. In the case of 
a cross section having double symmetry the four components of axial 
displacement are governed hy the following four equations: 






(80a) 









(80b) 









(80c) 



(80d) 


In the case of a cross section which is symmetrical about the 
x-axis the two components of displacement are governed by the following 
two equations: 






,2 

d w. 


dz^ 


% 




H 

GAg 


(8la) 



(8lb) 
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In the case of a cross section which is symmetrical about the y-axis 
the two components of displacement are governed hy the following two 
equations: 



V 

" 'GA^ 


(82a) 




d w 




dz^ 


5Ir, 


t/j 




T H 

GIc ’'tJ " GAg ’hJ 


(82b) 


Solution of Homogeneous Equation 

The solution of the nonhomogeneous problem may again be expressed 
in terms of the natural displacement modes which arise from the homo- 
geneous equation after a separation of the variables. The axial dis- 
placement nay be assumed as the product of a function of j by a 
function of z . 


wj = fjg(z) (83) 

For convenience the doubly symmetrical cross section will be considered 
with torsional loading (the antisymmetrical modes) . The homogeneous 
equation is obtained from equation (80b) by setting the right-hand side 
eq'ual to zero. 


+ f ^ ^ ^ = 0 (84) 

tJ 

If equation (83) is substituted into equation ( 84 ), the variables 
may be separated in the same manner as explained in section I. The 
spanwise and transverse distributions are found to be governed by the 
following two equations: 

= 0 (85a) 

dz^ 


^ ^ j - 0 


( 85 b) 
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2 

where = —5-. Unfortunately the method of solution for the stress 

modes which was used in section I is not particularly convenient here. 
The solution of equation (85h) may he obtained hy several different 
methods. Equation (85h) involves only one independent variable. In 
such cases, as has been previously noted, the use of generalized differ- 
ences offers no assistance in the solution of the problem although they 
do permit the writing of the equation in a condensed form. 

In order to solve equation ( 85 b) it is necessary to recognize that 
it represents a homogeneous system of linear algebraic equations which 
must be solved simultaneously. For simple cross sections having only a 
few stiffeners the eqmtions may be solved by direct algebraic methods. 
For sections having a large nimiber of stiffeners numerical methods of 
solution must be en^jloyed. In order to proceed to the solution of the 
original nonhomogeneous problem, it will be assumed that the natural 
stress modes have been con 5 >uted. 

The equations which govern the stress modes for other types of 
loading and symmetry conditions will now be stated briefly. For a beam 
with a doubly symmetrical cross section acted upon by horizontal loads, 
the stress modes are symmetrical about the x-axis and are governed by 
the following equation: 




X^a.f . 






(86a) 


For a beam with a doubly symmetrical cross section acted upon by verti- 
cal loads, the stress modes are symmetrical about the y-axis and are 
governed by the following equation: 


i(Wj) 


X^ajfj 


hA 




f. = 0 


(86b) 


If the cross section of the beam has single symmetry about the 
x-axis, the symmetrical modes are governed by equation (8&) while 
anti symmetrical modes are governed by the following equation: 









(87a) 


When the cross section of the beam has single symmetry about the y-axis, 
the symmetrical modes are governed by equation (86b) and the 
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antisymmetrical modes axe governed "by the following equation: 

If the cross section is completely unsynnnetrical the stress modes are 
governed hy the following equation: 

( 88 ) 


Solution of Nonhomogeneous Equation 

Before solving the nonhomogeneous equation it is necessary to prove 
the orthogonality of the stress modes and to introduce a law of normali- 
zation as in section I. Since the proof of orthogonality follows the 
same method used in section I it will he omitted. The condition of 
orthogonality and law of normalization nay he stated in one equation 
hy using a Kronecker delta. 




a.f . f . =8 

d jm Jn mn 


(89a) 


where 


8 ^ =-^ 


1 , 

0 , 


m = n 
m 5 ^ n 


(89b) 


In equation ( 89 a) fjjjj^ and fj^^ are the mth and nth modes, respectively. 
Equation (85h) may he written for the nth mode to obtain 




( 90 a) 


where 


"Tn ~ 




( 90 h) 
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Multiplying throng equation ( 90 a) by and summing over all 

stiffener areas give, after introducing equations (89), 

0 0 

( 91 ) 

Equations (89) and (91) may now be employed in solving tbe nonhomogeneous 
equation. 

A beam having a cross section with double symmetry and acted upon 
by torsional loading will again be considered. The axial displacement 
may be assumed to be given by a finite series as follows: 




(92) 


Substituting this formula into equation ( 80 b) gives, after interchanging 
the order of summation in the third term, 

y~ g A (3 A f ^ + I a f + 

^ n J \ t k Jn) G j 2 ^ jn Z_ Tn n Glc Tj 

Equation ( 93 ) must now be multiplied throu^ by fj^ and summed 
over all stiffeners. Reversing the order of summation in the first two 
terms gives 

^2 

5Z ^ 8n^‘jm(^J (PlAt^Jn) + | XI 5Z ® j^Jm^jn + 

n J ' n J dz 


^Tm P T P 

Ic A- 

Introducing equations (89) and reananging terms gives 


E 

G . 2 
dz 


H 

n 


Sn 






V 

GI 


Substituting equation (91) and changing the subscripts m to n give 
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or 


G dz2 ■ ^ ^ GIc 

2 _ 
gn “ei 



(9M 


Equation (9^)^ ■with the end conditions, gi'ves the spanwlse distribu- 
tion of each stress mode and thus coD 5 )letes the solution for the 
displacements . 


Similar equations may he developed for the spanwise distribution 
of stress modes ■with other types of loading and symmetry conditions. 
Since the mathematical method of development of the equations is the 
same in all cases, the equations ■will he merely s'tated -without proof. 
For a doubly symmetrical section siibjected to horizontal loading the 
stress modes are det<.rmined by 

^^gn 2 

where 

^Hn = ^^^c) 



For a doubly symmetrical section subjected to vertical loading the 
stress modes are determined by 




dz 


n 2 
2 “ hi gn 




(94d) 


where 


"Vn 




’'Vj^jn 


(9he) 


When a cross section has single symmetry about the x-axis, the 
symmetrical stress modes are governed by equation (9^t)) while the 
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antlsymmetrical inodes are governed by the equation 


is ^ W W 

dz2 ~ "EI^, - EAv 


( 9 ^f) 


If tbe section has single symmeti^ about the y-axis, the symmetrical 
modes are governed by equation ( 9 ^) 'wbile tbe anti symmetrical modes 
are governed by tbe equation 


^ P V ^ 

^^2 Sn - "EIc ' EAh 


(9^) 


¥ben tbe cross section is con^jletely unsymmetrical, tbe stress modes 
are governed by tbe following equation: 

- W V V 

d^2 ^ EAh eav 

A convenient method for calculating tbe section constants 
and Cyj^ is given in appendix D. 



Evaluation of Difference Term 

In eacb of tbe equations for transverse stress, or displacement, 
modes there appears a difi^rence term of second order. This term must 
be expressed in algebre^ic form before tbe stress modes can be determined. 
Tbe difference term has the form Aj(PT5j!!^f j]. In converting this dif- 
ference quantity to algebraic form it is necessary to express all sub- 
scripts in terms of stiffener numbers. Tbe subscripts referring to 
Vebs must be eliminated. 

In figure 4 (b) there is shown a segment of an idealized stiffened 
section in which web a and web b are connected to the jtb stiffener. 

Tbe first-order difference quantity A]j^f j is expressed, over web a, 
in tbe following manner: 

- "3-1 

Tbe quantity Pij^fj becomes, for web a. 
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For web b this quantity becomes 


In tbe above equations it should be noted that the single subscript on 
the relative warping stiffness 3 has been replaced by a double sub- 
script corresponding to the stiffeners at either end of the i>artictilar 
web being considered. 

The second-order difference quantity may now be expressed in the 
following algebraic form: 


where 


The quantity EjP is the sum of the relative warping stiffnesses of 
the webs which connect to the jth stiffener. 

It may be of some interest to note the algebraic form of the dif- 
ference term when more than two webs connect to the stiffener as may 
occur in multicell sections. In figure 6 there are three webs connected 
to the Jth stiffener. The assumed positive directions for all flow- 
type quantities are indicated on the webs. The difference term, for 
this case, becomes 

= Pjm^m + Pj(j+l)^j+l - + Pj(j-l)^j-l (96a) 

where 

¥ = ^96b) 


Modified Stiffnesses for Symmetrical Cross Sections 

In any idealized stiffened cross section the number of stress modes 
will be equal to the number of stiffeners. For a section with n stif- 
feners there will be three primary planar modes and n - 3 secondary 
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self-equilitrating modes. For an tinsymmetrical section equation ( 88 ) 
corresponds to n linear algebraic equations wMcb may be solved to 
obtain tbe stress modes. If a cross section has an axis of symmetry^ 
the symmetrical and anti symmetrical modes may be conq)uted from two 
independent smaller systems of equations. Tbe conditions of symmetry 
may be conveniently incorporated into tbe equations by using modified 
stiffness values for webs or stiffeners which lie on an axis of 
symmetry. 

As an exanq)le a section may be considered which is symmetrical 
about the y-axis and which has a web on the y-axis as shown in figure 8 . 
The symmetrical modes are governed by equation ( 86 b) which may be 
written for every stiffener in the cross section. However, it is only 
necessary to write this equation for the stiffeners on one side of the 
axis of symmetry. For a symmetrical mode the axial displacements will 
be uni form over the central web as shown in figure 8 (a) . The axial dis- 
placement and normal stress will be the same in stiffeners 3 and k. 

If equation ( 86 b) is written for stiffener 3t and equation (95a) is 
substituted, the following equation may be obtained: 

7 

-P3l|f4 + (233)^3 - ^32^2 - 

Assuming that a system of equations wooild be formed for all stiffeners 
on the right side of the axis of symmetry, it becomes essential to 
eliminate fj^. from equation (97)* Since f}^ is equal to f^ it may 

be replaced by f^ to obtain 

J 

^32^3 " ^ 32^2 “ ^ ®'3^3 ( 98 a) 


since 


+ ^32 (98b) 

Equation ( 98 a) could have been written immediately by omitting the tern 
containing the stress value on the left side of the axis and assuming 
the warping stiffness of web 3 ^ to be zero as shown on the idght-hand 
side of figure 8 (a) . 

In a similar manner the anti symmetrical mode may be considered as 
illxistrated in figure 8 (b) . In this case it is foimd that the standard 
stiffness value must be multiplied by a factor of 2 as indicated on the 
right-hand side of figure 8 (b) . 
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The cross sections may have a stiffener on the axis of symmetry as 
shown in figure 9 * The equation which is associated with this stiffener 
may he formed, for the symmetrical mode, hy multiplying the stiffener 
area hy one-half as shown in figure 9 (a) and omitting the term con- 
taining the stress in the stiffener on the left side of the axis. If 
there is a geometric discontinuity at this stiffener a factor of one- 
half must he applied to this quantity. If there is an internal weh 
lying on the axis its stiffness must he multiplied hy one-half. 

Although the use of modified stiffnesses does provide some convenience, 
the advantage of their use is not large since the symmetry conditions 
are easy to apply. For a section with double symmetry it is only neces- 
sary to consider the naterial of the cross section that is contained 
within the first quadrant. 


Algebraic Equations for Stress Modes in a Doubly Symmetrical Section 

The bending and torsion modes for a doubly symmetrical section are 
governed hy equations (85b), (86a), and (86h). An equation for the 
doubly symmetrical modes, which are due to axial force, may be obtained 
from equation (80a) by a separation of the variables. An algebraic 
form of the equations for the stress modes may be obtained by replacing 
the second-order difference tenn by its algebraic eqiiivalent as given 
by equations ( 95 ) or (96) . 

For puiposes of illTistration a single-cell section may be con- 
sidered in which each stiffener will be connected to two webs. Typical 
equations for the Jth stiffener in each stress mode appear as follows: 


“Pj(j+l)^j+l ■ ^^(J-D^J-I + ^ = 0 (99b) 

y 

"Pj(j+l)^J+i ^ ^ ° ( 99 d) 
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Each of the above equations is a typical equation of a system of equa- 
tions governing a particular type of mode. The number of equations in 
each system depends upon the number of stiffeners in the first quadrant 
having finite values of stress. The process of formation of similar 
equations for sections having single symmetry or without symmetry is 
obvious. Examples of the application of these equations to specific 
cross sections are given in a subsequent section. 


Matrix Forms of Eqtiations for Stress Modes 


The stress modes for an unsymmetrical section are governed by 
equation (88) . In order to express the corresponding system of algebraic 
equations it is convenient to introduce the notation of matrix algebra. 

In order to make clear the meaning of each matrix symbol it is desirable 
to show the matrices in expanded form. In order to avoid the necessity 
of writing out large natrices the expanded matrices will be written for 
an unsymmetrical four-flange section. The contracted matrix equations, 
however, will have conqjlete generality. 

Before writing the equations it is useful to note that the summa- 
tions which occur in equation (88) may be expressed as scalar products 
of vectors. If [f ] indicates a column matrix, or colinnn vector, 
then [f ]’ will indicate a row matrix or row vector. The prime is used 
to indicate the transpose of a mato. If [f] is a vector of stress- 
mode ordinates, then [f ] and [f] for a four- flange section are 
given in expanded form as 



[f]' . f2 fj 


( 100 ) 


Vectors of the geometric discontinuities associated with the flanges 
are defined by 



(!■ 


^T1 ^T2 ’'tS 



fll ^H2 ’'h3 
I^Vl ^V2 ’'V3 


(101a) 

(101b) 


(101c) 
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The summations -which are con-talned In equation (88) may now he expressed 
hy the following scalar products: 


= H'H 


(I02a) 


^ = f'v] 0 

Equation (88) my now he -written In the following form: 


(I02h) 


(I02c) 






= 0 (X03) 


Substituting equation (95a) and transposing the term containing 
gives 




In this equation the coefficients of the first three terms are 
warping stiffnesses- These coefficients, for a fo-ur-flange heam, my 
he arranged to form the follo-wlng mtrlx: 




-Pl2 

-P21 


0 

-P32 

-\i 

0 


0 

~Hh 

"^23 

0 

Z3P 

-334 

-^3 



(105) 
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The rows of this matrix contain the warping stiffnesses as they would 
appear in the four algebraic equations for the stress modes. It may 
also be noted that the subscripts correspond to the position of the 
element in the matrix. This will be true for any cross section. The 
diagonal element of the jth row is the sum of the stiffnesses of all 
webs which connect to the jth stiffener. The nondiagonal nonzero 
elements of the jth row are the individml stiffnesses of the webs con- 
necting to the jth stiffener. From these considerations it is seen 
that this matrix -can be easily formed for any cross section from the 
web dimensions. 

On the right-hand side of eqmtion (l04) there appears the stif- 
fener area. These areas may be formed into the following diagonal 
matrix: 


a]^ 0 0 0 



0 0 
0 a^ 0 


0 0 0 



(106) 


Using the definitions of equations (105) and (106), the system of equa- 
tions, corresponding to equation (l04), maybe written in the following 
matrix form: 


B f 


HH'H - 4HHT'] - ■ >'[*]['] 

(107) 


When a column vector is premultiplied by a row vector the resxilt is 
a scalar number. Hence the product may be said to be a scalar multipli- 
cation. Exan^jles of such products are given in equations (l02). In 
equation (lO?) it is seen that there are three terms in which column 
vectors of geometric discontinuities are postraultiplied by row vectors 
of the same geometric quantities. When a column vector is postraultiplied 
by a row vector the result is a square matrix and the product may be 
called a matrix product of two vectors. It is convenient to introduce 
matrix symbols for these products as follows: 



(108a) 

(108b) 
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(108c) 


In expanded form these matrices appear as follows for a four-flange 
■beam: 





2 

7^1 

^T 1 ^T 2 

^Tl^TB 

^Tl^Tlt 

^T 2 ^T 1 

7 2 

7q,2 

^T 2 ^T 3 

7^2^T4 


^T 3 ^T 2 

7 2 
75,3 

^tb’'t4 

y y 
Tl^ T1 

y y 
T4 T2 

y y 
t4 TB 

1 

f7 2 

^H 1 ^H 2 

%L^H3 

'^m7Ek 

^H 2 ’'h 1 

7 2 

^H2^H3 


^H 3 ^E 1 

’"hb’'h 2 

2 

’hb 

y y 
H 3 h4 

Whi 

Wh2 


Y 2 

\k 

CVJ 

^V 1 ^V 2 

'^yJy3 

’Vi’'v4 

7y2^Y1 

ro 

ro 

7y2’Vb 



^V 3 ^V 2 

7 ^ 
V3 

7^7„, 
V3 Vif 



WVB 

7 2 


(109a) 


(109b) 


(109c) 


The formation of these nstrices is obvious hy inspection since the 
subscripts correspond to the position of the element in the matrix 
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Equation (107) may now "be written in the following form: 



The vector [ f ] may he factored out on the left-hand side to obtain 

' (H - 0 

V. ^ XJ. ' J 

The matrices within the braces may be combined to give this equation in 
the form 

[C][f] = X^A][f] ( 110 a) 

where 



( 110 b) 


The matidces [A] and [C] are dependent solely upon the geometric 
properties of the cross section. Convenient methods for computing the 
elements of (T^ shown in appendix C. The 

only infonnation that is required for the conqnitation of these latter 
matrices is the coordinates of the stiffeners with respect to the prin- 
cipal shear axis. Thtis the properties of the cross section which enter 
into the matrices [C] and [A] are the areas of the stiffeners, the 
coordinates of the stiffeners, and the waiping stiffhesses of the webs. 

A few comments concerning numerical methods of solution of equa- 
tions (110) may be worth while. It would seem, by analogy with vibration 
problems, that the method of matrix iteration coiild be applied after 

multiplying through the equation by [A] . However, unfprtxmately, the 

process proves to be divergent. If the matrices are reduced to eliminate 
the planar modes, the iteration process will continue to be divergent. 
There are several numerical methods which may be used to solve these 
equations. However, a study of the relativfe merits of such methods is 
an extensive project in itself ahd will not be treated herein. 

Equations governing the stress modes in sections having single or 
double symmetry will now be considered. From arguments similar to those 
used to derive equations (llO), two matrix equations maybe obtained 
for the symmetrical and antisymmetrical modes for a cross section which 
is symmetrical about the x-axis as follows: 
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vliere 




( Symmetrical) 
(Anti symmetrical) 



(111a) 

(lilt) 

(111c) 

(llld) 


If tlie cross section is symmetrical about tlie y-axis, eqxiations (ilia) 
and (lllb) may again be used witb the following definitions of the 
coefficient matrices : 


EG ' [?] - ^Ev] (112a) 

EG = H - ^EG - ^EG 

If the cross section has dotible symmetry, there are four types of modes 


as detennined by the following 
where 

M- 

D= J - [B] 

K * H 


matrix equations; 


(Symmetrical) 

(113a) 

(Anti symmetrical) 

(I13t) 

(Symmetrical about x-axis) 

(113c) 

(Symmetrical about y-axis) 

(113a) 



(I13e) 

1 

" Ic 

EG 

(I13f) 

Att 

Eh] 

(I13g) 
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ty] = [®] - ai3n) 

For simple cross sections algetraic solutions of equations (ill), (ll2), 
and ( 113 ) inay readily derived. A few exanq)les are given in the 
following section. 


Examples of Stress Modes and Section Constants 

A few simple exanq)les will be solved for the stress modes and 
section constants. Fotir rectangular sections, having various numhers 
of stiffeners, and also a trapezoidal section, will he considered. The 
rectangular sections will have doiible symmetry while the trapezoidal 
section will have single symmetry. 


Fovir-flange section with double symmetry .- The dimensions • of the 
section are indicated in. figure 10(a) and the web stiffnesses are shown 
in figure 10(b) . Since the section is doubly symmetrical it is only 
necessary to consider the first quadrant with modified stiffnesses as 
shown in figure 10(c) . In the case of a four-flange section there are 
three planar modes and one secondary mode. The secondary mode is gov- 
erned by equation (99b) which takes the following form: 


L^)h - 




Tl^l) 




The coefficient of the first term is given by 

The geometric discontinioity is given by (see appendix D) 

7ti = 

‘ = t(^^ - \) 

The central moment of inertia is given by (see appendix C) 


(lll^) 


(115) 


(116) 






(Pb + Ph) 


(117) 
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Inspection of equation (ll^) stiows that the stress is arbi- 

trary. Its value must he determined hy the normalizing condition as 
e:jq)ressed hy equations ( 89 ) • The following formula is obtained; 



( 118 ) 


The remaining stresses are obvious from the conditions of symmetry. The 
only remaining unknown in equation (ll4) is the characteristic numiber X. 
Siibstituting equations (115), (H6), and (II 7 ) into equation ^l4) leads 
to the following formula for X^: 


,2 

This formula was originally given hy Grzedzielski (reference 2). 


(119) 


In order to determine the spanwise variation of the stress mode it 
is necessary to detennine the section constant Cfji. This constant is 
given hy the following formula: 






( 120 ) 


It is also of interest to note the following section properties which 
are obvious in this case: 


% = 2h1>b 
Ay = 2htj^ 


( 121 ) 


Six- flange section with double symmetry .- This section is illus- 
trated in figure 11. There are six stress modes of which three are 
primary and three are secondary. The doubly symmetrical modes will he 
first considered. Only the flxst quadrant, as shown in figure 11(b) , 
need be considered with modified stiffnesses as shown. Corresponding 
to equation (99a) two equations may be written as follows: 

(^12 " ®’ l ^ )^1 ” ^ 12^2 “ ® 


(I22a) 
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-^12^1 + (Pl2 - = 0 (122b) 

Tbe determinant of tbe coefficients must vanish. 

(Pl2 - ^ 1 ^^) 

-P12 

Since the primary planar mode has the characteristic value X = 0, 
the expansion of the above determinant must lead to an equation -which 
is linear in \r. This equation may be immediately solved to give 




12 




= 0 


(123) 



(124) 


The determination of the secondary stress mode consists merely of 
sol-vlng for f^ in terms of f^. This may be done by substituting 

equation (124) into equation (l22a) or (l22b) . It may also be done by 
using the orthogonality condition bet-ween the primary and secondary 
modes. After normalization the stresses thus obtained are 



( 125 a) 



(125b) 


Four stress modes renffiiln to be determined. There -will be one 
doubly anti symmetrical mode. In this mode the stresses in stiffeners 2 
and 5 will be zero. The section -will thus behave as a four-flange sec- 
tion under torsional loading and the formulas which have been previously 
gi-ven may be -used to determine this mode. If the beam is subjected to 
horizontal bending there -will be one planar mode which may be computed 
from the flexure formula. 

Under the action of vertical bending t-wo stress modes -vd.11 arise. 
These modes -will be symmetrical about the y-axis and anti symmetrical 
about the x-axis. The essential stiffnesses are shown in figure 11(c). 
The stresses are governed by the following two equations: 
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- Pl2^2 





(126a) 


“ Pl 2^1 ■*■ Pl 2^2 ~ ^ "^.^2 (126I3) 

l-Hieii ttiese eguations are solved the normalized values of f]_ and fg 
for the secondary mode are fomd to he the same as given hy egua- 
tions (125). The value of for the secondary mode is given hy 

equation ( 124 ). In order to calculate the spanwise variation of the 
secondary mode it is necessary to determine the section constant 

Using equation ( 94 e) this constant is foimd to have the value 


Cy2 = - 



(127) 


Six- flange sections vlth double symmetry and two cells .- This cross 
section is Illustrated in figure 12 . It differs from the previous cross 
section hy having an internal vertical shear weh. It may he readily 
determined from symmetry conditions that five of the six stress modes 
for this section are identically the same as for the previous section. 
The mode which is different is the secondary mode due to vertical 
bending. Since the mathematical methods of solution are the same as 
before, the final formulas will he stated withotrb derivation for the 
secondary vertical bending mode only. 


The characteristic number is given hy 





(128) 


The stiffness quantities Z ]_3 and Z2P ^.re modified stiffness values 
as shown in figure 12 (c) . They are given hy the formulas 
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The normalized formulas for the stresses are again given hy equa- 
tions (125) . The section constant Cy2 given hy 


Cy2 =- 


^s 




(129) 


Eight-flange section with double symmetry .- This cross section is 
illustrated in figure 13(a). There are two modes having each type of 
symmetry. The stresses and characteristic number for the secondary mode 
having double symmetry are obtained from equations (l24) and (125) by 

^2 

substituting a£ for The two doubly antisymmetrical modes are 

governed by the following equations: 

y 

(2lP)fl - = VaA (UOa) 


■^ 12^1 (^ 2^)^2 “ ^%^2 ^^ 30b ) 

The characteristic equation is as follows: 



(131) 


Both modes are secondary and it is necessary to solve a quadratic 
equation in Consequently no simple formiilas can be giyen for 

these two modes. 


In the case of vertical bending there is one primary mode and one 
secondary mode. For this secondary mode the stresses and characteristic 
numbers are agalng^obtained from equations (124) and (l25) by substi- 
tuting &2 ~2’ associated section constant Cy2 is given by 
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Cy2 - - 


^•bl8 1 ^ 


(132) 


The equations vhich govern the horizontal tending modes are as follows: 

(SiP)fi - ^ ^1^1 (133a) 

J 


“^ 12^1 (^ 2 ^) ^2 = 

2 

The characteristic equation is linear in X and leads immediately to 
the following fonmila: 


X 


2 




(m 


Orthogonality between the primary and secondary mode is expressed by 
the following equation: 


= ° 


or 



If the above formula is substituted into the normalizing equation, it 
is found that the normalized stress formulas are as follows: 



( 135 a) 


( 135 ^) 
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The associated section constant is given hy 


Cg-2 - 


4t^pXp 

\f^ 



( 136 ) 


Trapezoidal section .- This section is assumed to have single sym- 
metry as shown in figure l4. The section properties % and Ay may 

he immediately determined to he as follows: 


Ah = 2tj2Li2 cos^ai2 


(137a) 


Ay = + ^23-^3 sin^otQj2 (l 3 Tt>) 

The remaining section properties cannot he determined until the prin- 
cipal shear sixes are located. Because of symmetry it is only necessary 
to locate the origin of the axes. This is determined from equa- 
tion (71b). This equation provides the following formula for the 
coordinate of the first flange: 

V " :^(‘^23^3 ■ ■*^12^lJ+ “is) 

This formula may he obtained immediately hy direct substitution into the 
formulas given hjr Kirste (reference 5)> Drymael (reference 8 ), or Duherg 
(reference 18 ). The rsidii to the wall segments and the geometric dis- 
continuities may now he computed. 

Two of the stress modes are symmetrical about the x-axis. These 
are primary modes due to axial force and horizontal holding. The 
remaining two modes are antisyxmnetrical about the x-axis. One of these 
is a primary mode and the other a secondary mode. The antisymmetrical 
modes are governed hy the following equations: 

(V)^l " ^12^2 " ’'t/o - ^ ^1 (139a) 


- Pl 2^1 + (^^)^2 - ^ ^ ^ t/j - - a 2^^^2 ( l 39 b ) 
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The characteristic equation is linear in X a3od yields the following 
formula: 




27ti^ _ ^ a 4 3 _ 


a^V 2- 


(1^) 




In the ahove equation the summations ^2^ contain modified 

stiffnesses as illustrated in figure llf(h). 

The stress values are most easily obtained from the orthogonality 
condition vhich is expressed as follows: 


a^y^fl + ^ 


or 



( 1 ^ 1 ) 


The normalizing condition then gives the following specific stress 
formulas : 



(ll)-2a) 


(l42b) 


The section properties 


and are given by the formulas 



7viy2 







(ll^3b) 
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vhere 


7vi = ti 2 sin 

(ll^3c) 

7^2 = ~^23 ~ ^12 “12 


T1 ^ ^ 12^12 " ^ll<-^lll- 

(li^3e) 

T2 "" ^23'*^23 “ ^12^12 

(l43f) 

coHCLUDnra eemarks 



A method of solution has heen obtained for dete rminin g secondary 
stresses in thin-walled beams of arhitrary cross section. Monocoque 
sections with comers require additional conditions of continuity at 
the comers and have not heen considered. The introduction of general- 
ized difference equations for the semimonocoque sections appears to 
offer definite advantages and may prove to he useful in the analysis of 
other types of stmctures. By using these difference equations the 
analysis for stiffened heama can he carried out in a manner exactly 
parallel to the method of analytic functions for a single-cell 
unstiffened heam. 

Future research study should consider elastic rihs for wings, 
elastic stiffening rings fbr noncircular fuselages and the thick- 
walled wing without rihs. It is possible that the natural stress modes 
will he useful in such prohlems. An extension of the present theory to 
swept beams is also a currently important problem. Although of lesser 
practical Importance, the effect of secondary stresses on vibration 
frequencies may now he determined. Practical numerical methods of 
computing the natural stress modes of stiffened sections should he 
developed. A comparison of all available test data with accurate, as 
well as approximate, solutions for secondary stresses would he of 
considerable Interest. 


California Institute of Technology 

Pasadena, Calif., January 21, 1951 
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APPENDIX A 

DEFIiKTIONS DUE TO PLANAR STRESS DISTRIBUTION 


Equations ( 23 ) give the following relations between stresses and 
deflections : 


d u _ ^ 
dz^ 

d^ _ Pv 




cos a ^ ds 


Scr 

sin a 


ds 


(Al) 


(A2) 


The flexure formula defines a planar distribution of stresses as follows 


a 


p MyX M^y 

"A;."' ly'" Ix 


(A3) 


It is necessary to express the centroidal coordinates x and y 
in equation (A3) in terms of the coordinates x aud y with respect 
to the principal shear axes. The relationship between the two coordi- 
nate axes is shown in figure I 5 . Ercm the geometry of this figure it 
may he shown that the coordinates of a point are related by the 
following equations: 


X = X + X cos 0 - y sin 0 
o 


Uh) 


y = y + X sin 0 + y cos 0 
•' •'o 


(A5) 


The differentials dx, dy, and ds are related to the angle a 
by the following formulas: 


dy . ' 

^ = sin a 
ds 


dx 

ds 


= cos a 


(a6) 




It is necessary to determine ^a/be to substitute into equa- 
tions (A1) and (A2). Differentiating equation (A3) gives 
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^ ^ ^ ^ 

“I ds ^ I ds 

y X 

Differentiating equations (a 4-) and (A5) gives 

■^ = T- cos 0-3^ sin 9 
ds ds ds 


(A7) 


(a8) 


dy dx dy 

di" = dJ ® + di’ ® 


Substituting equations (a 6) into eqmtions (A8) and CA9) gives 


dx 


3 — = cos a cos 9 - sin a sin 9 
ds 


. (A9) 


(AlO) 


ds 


= cos a sin 9 + sin a cos 9 


(All) 


Substituting equations (AlO) and (All) into equation (A7) gives 


5cr „ \ Mx/- \ 

^ = ^Mcos a cos 9 - sin a sin 6) + s^cos a sin 9 + sin a cos 9) 

(A12) 

The integral contained in equation (Al) may now be evaluated as follows: 
cos a ^ ds = cos a(cos a cos 0 - sin a sin 0) ds + 


hr 

Mx 


^wt^ cos <x(cos a sin 0 + sin a cos 0) ds 


(AI 3 ) 


Two of the four integrals on the right-hand side vanish because of the 
definition of the location of the principal shear axes as expressed by 
equation (l5c). The integral takes the following value: 




ba ^ ^ 

cos a ^ ds = Ag cos 0 + rjr- Ag- sin 0 

os ly • 


(Alh) 
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Similarly tlie integral in equation (A2) may te shown to have the 
following value: 



sin a 


ScT 


ds 



% 

sin 0 + =— Ay 
rx 


cos 9 


(AI5) 


Equations (Allf) and (AI 5 ) may now he substituted into eq-ua- 
tions (A1) and (A2) to obtain the following formulas for deflections: 


d^u PH My. cos e Mjj. sin 9 
^ = -G% - -E^i; 


(AI 6 ) 


Py Sin 9 cos 9 

^ = E^ 


(A 17 ) 


These equations are given in the text as equations (25). 
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APPENDIX B 

efeex:ts of symmetry of cross section 


The axial displacements in an tmsymmetrical cross section are 
governed "by equation ( 20 a). From considerations of symmetry it is 
possible to replace this equation by four equations of simpler form. 

It is well-known that solutions to structural problems involving one 
independent variable can be divided into symmetrical and antisymmetrical 
parts when the structure is symmetrical. In the case of a beam with a 
cross section having double symmetry, since the transverse distribution 
of stresses or displacements on a cross section depends on two inde- 
pendent variables , the solutions may be divided into four parts or com- 
ponents. The four con5)onents of a function will be indicated, in the 
present appendix only, by subscripts defined as follows: 


a antisymmetrical about both axes 

s symmetrical about both axes 

X symmetrical, about x-axis 

y symmetrical about y-axis 

The stresses and displacements which occur on a cross section may 
be expressed in terms of their components as follows: 


= ^a 

+ 


+ 

•^x 

+ 


Q 

II 

Q 

+ 

0 

+ 

0 

+ 

0 

a 


S 


X 


y 

% = \ 

+ 


+ 


+ 


a 


a 


X 


y 

w = Wa 

+ 

^s 

+ 

^x 

+ 

^y 


J 


(Bl) 


Three geometrical quantities r, cos a, and sin a, enter into the 
analysis. For a section with double symmetry these quantities have the 
following symmetry properties : 


r 


Symmetrical' about both axes 


cos a 


Symmetrical about y-axis 


sin a 


Symmetrical about x-axls 
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The symmetry properties of the derivatives of these geometric quantities 
are readily seen to he as follows: 

dr/ds > Antisymmetrlcal about both axes 

d cos a/ds > Symmetrical about x-axis 

d sin a/ds > Symmetrical about y-axis 


It is now necessary to determine the symmetry properties of the 
second derivative of w with respect to s. The four components of 
the second derivative may be indicated as follows: 


Ss^ 



(B2) 


If a function having double symmetry is differentiated once, a function 
having double antisymmetry is^ produced, and vice versa. If a function 
which is symmetrical about the x-axis is differentiated once, a function 
is obtained which is symmetrical about the y-axis, and vice versa. 
However, if a function has a particular type of symmetry, the second 
derivative of the function will have the same type of symmetry. Con- 

a^w 

sequent ly the four components of — ^ can be written in terms of the 
four components of w as follows: 






(B3) 


In addition it is necessary to consider the integrals which appear 
in equation C 20 a) . After an integration by parts the integrand of each 
integral contains the product of a geometric property and the nxi .al 
displacement. The only component of the axial displacement which con- 
tributes a finite quantity to the value of a given integral is that 
cdnq)onent which has the same symmetry as the geometric property 
appeairing in the integral. From this consideration the three integrals 
are seen to have the following values: 
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(Bl^) 


It is now possible to consider tbe symmetry properties of each 
term in equation (20a). The first two terms have four cong)onents as 
defined by equations (Bl), (B2), and (B3). Each of the remaining terms 
h»n only one type of symmetry which may be determined from the symmetry 
properties of the derivatives of the geometric properties of the cross 
section. The sum of the terms on the left side of equation (20a) 
having a particular type of symmetry must be equal to the term on the 
right side having that seime type of symmetry. This permits the writing 
of four equations corresponding to the four types of symmetry as shown 
below. 


Symmetrical about both axes: 


^2 ^2 

o ws E» o Ws 


" “ r> 

__ . — s= U 




Antisymmetrical about both axes: 

T? 1 


Symmetrical about x-axis 


E» ^^a 1 dr^^ ^a ^ T dr 

G" ^;2“ ■ J~ dBj ^ 'ST ^ ~ "GIp ds 


(B5) 


(b6) 


^^x E ' ^^x 

2 G ^ 2 


1 d cos a 
Aj SB 


3b 3z 

Symmetrical about y-axis: 


a-^ds =- 


H d cos a 




ds 




z + ® 


Ss' 


1 ^Blna^ Bina^as 


V d sin a 


ds 


(B7) 


(B8) 
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A corresponding set of four equations for normal stresses is readily 
ottalnatle. 

The displacement dlstrihution governed, hy equation (B5) is due to 
the action of an axial force only. Methods for the solution of this 
elliptic equation are well-known. The problem is mathematically equiva- 
lent to a flat rectangular plate with transverse stiffeners of infinite 
stiffness and with sidesway prevented.. Equations (b 6), (B7), and (b 8) 
are given in the text as equations (2?) with the subscripts omitted, for 
convenience. The sTibscripts which appear in equations (B5) to (b 8) 
show clearly the particular component of the displacement which is 
associated with each type of external loading. Erom the differential 
equations which relate cr, and 5 to w it is possible to deter- 

mine which component of t, a, or I occurs with a particular type of 
loading on a beam with doubly symmetrical cross section. 

Men a cross section is symmetrical about one axis only, the solu- 
tion for an unknown may be divided into symmetrical and antisymmetrical 
parts only. For a cross section which is symmetrical about the x-axis 
only, the geometrical properties of the cross section have the following 
symmetry properties: 

r > Symmetrical aboiit x-axls 

cos a > Antisymmetrical about x-axls 

sin a > Symmetrical about x-axis 

dr/ds > Antisymmetrical about x-axis 

d cos a/ds > Symmetrical about x-axis 

d sin a/ds ^Antisymmetrical about x-axls 

By using arguments similar to those used for the doubly symmetrical 
section, it is found that equation (20a), for the general cross section, 
separates into the following two equations for sections which are sym- 
metrical about the x-axls: 


>2 
o w. 


sx 


8s^ 


E' 

G 


;^2 
o w. 


sx 




1 d cos a 


ds 




COB a 


8wsx 

ds 


ds = - 


H d cos a 




ds 


(B9) 
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ax E* 
2 G 




ax 




1 y. ^ax 

I_ ds y ds 


ds 


1 d sin ay^j. 
Ay ds ™ 


^ax 

sin a — r — ds = 
os 


T dr V d sin a 
GZ^ dB CcA-y cLs 


(BIO) 


In these equations the suhscripts sx and ax indicate symmetry and 
antisymmetry, respectively, about the x-axls. Equations (B9) and (BIO) 
are given without subsciipts as equations (29) in the text. 

In a similar manner it may he shown that, for the case of a section 
which is symmetrical about the y-axis, the general equation separates 
into the following two equations: 


^^sy E ' ^^sy 



d sin a 
ds 



sin a 




sy 


6b 


ds 


V d sin g 
GAy ds 


(BU) 


^2 ^2 
° ^ay ^ ^ ^ay 
as2 G ^^2 




1 d cos a 
Ag ds 



COB a 



ds = 


T dr H d cos a 
GI^ ds ~ GAg ds 


(B12) 


These equations are given without subscripts in the text as equations (30) • 
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APPENDIX C 

SHEARING SKITION PROPERTIES FOR STIFFENED SECTIONS 


Referring to figure l6(a), the kth veh is shown to extend from 
stiffener a to stiffener h. The length of this weh 1^ heis a hori- 
zontal componoit given hy Xg^ - x^, or and a vertical com- 
ponent - y^, or trigonometric functions of the 

angle ojj. may thus he written as 


■ cos 




(Cl) 


= lt 


(C2) 


The shear resistant areas are defined hy equations (72b) and (72c). 
They may thus he computed from the following formulas: 


% = 51 = 5Z (C3) 

k k ' ' 


Ay = sln^o^j. = (Afcyj)Sk (C 4 ) 

k k 

In order to develop a convenient formula for computing the central 
moment of inertia it is necessary to define the "central area" which 
may be associated with the kth web. The central area is Illustrated 

in figure l6(b) . It may be defined an the area of the triangle formed 
by the kth web and the two radial lines drawn to the origin of the 
principal shear axes from the ends of the kth web. The raditis is 

related to the coordinates of the stiffeners by the following equations: 


rij. = Xg sin - yg cos Oj^, 
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SubstltTiting equations (Cl) and (C2) gives 


2^ ya 


Lk 


- |("1> - 


' - ^a^b) 


(C5) 


The area Afc may nov he related to the coordinates of the stif- 
feners hy the following equations: 

= ^k\ - ^a=^b - 

This formula may he written as a determinant. 


2Av = 


yb 


(C?) 


The formula for the central moment of inertia is given hy equation (72a) 
and may now he cojnputed from the following formula: 


ic - H = H {^)\ 


(c8) 
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APPEHDIX D 

GEOMETRIC DISCOHTINUITIES OF CROSS SECTION 


Reference may be made to figure ^(b) where the Jth stiffener lies 
between web a and web b. The geometric discontinuites, which are 
defined at the stiffeners, are given by equations (77)* The discon- 
tinuity 7„. is defined by 

’’ej - >\) 

This may be written as 

= t^ cos - t^ cos 
Substituting equation (Cl) gives 

?HJ = - ^ 1 ) - - =<3-]) 

From equation (D2) it is seen that the column vector [7^] as defined 
by equation (lOlb) may be computed from the following matrix formula: 


M = ">3) 

where 

f 

H = 1^1-2 

The matrix [b] is defined for a four-flange section in equation (105). 

Equation (D3) is also applicable to multicell sections. From similar 
arguments the discontinuity 7^^ is fotmd to be given by the following 

formula: 


(Dh) 
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The column vector 



is given hy the equation 


where 



(D5) 


H = 


The discontinuity 7^^ is defined hy the following formula; 


=Ap_ = A/tr^ 
'Tj j'^Tk k kJ 

Introducing the central area Ajj. gives 


(d6) 




= A 


|("Pk\) 


(D7) 


The square matrices [T^ most readily computed 

from their definitions an given hy equations (108) after the column 


vectors 


r “1 


T 


7^, and 7^^ have been computed. 


The section constants ^Hn^ *^Vn^ corresponding to the 

nth stress mode, may he immediately computed from the formulas 


^Tn “ ^ ’^Tj V 

^Hn " ^ 


Vn 


’'Vj^on 
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APPENDIX E 

SUMMATION BY PARTS FOR gTIFFEEED SECTIONS 


The summation hy parts which is to he justified is given hy equa- 
tions (78) and hsLS the following form: 


The proof will he given for a very simple multicell section. The sec- 
tion to he considered is shown in figure 17. It has two cells and four 
stiffeners. The wehs and stiffeners are numbered and the positive 
direction of flow-type quantities is shown. 

The first summation of equation (El) may he expanded in the 
following manner; 

^ ' PiVj + P1.V3 * ‘ 

Pl("l - + P2(W2 - vj + P3(Wj - Vg) + - W3) + P5(»3^ - Wj) 

The terms on the right- hand side may nov he arranged in the following 
manner: 


^ pA-'j ' - Pa + P5) 

The summation which occurs on the ri^t-hand side of equation (El) 
may he expanded as follows: 


+ Vg(p2 - P3) + 

pj (12) 
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= Vj^(p2 - Pi - f^) + ''2(P3 




(E3) 


If equation (E 2 ) is now con5>aTed with equation (E 3 ), it is seen that 
equation (El) is valid. It is clear that the argument could he applied 
to sections having any finite number of cells or stiffeners since no 
new element would he hrou^t into the problem hy increasing the number 
of cells or stiffeners. 
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Figure 1. " Coordinate system . 
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(a) Dif-feren+ial elernent. 



(b) Axial forces. 



(c) Tangential forces- (d) Normal forcee . 

Figure 2.- Forces acting on a wall elernent. 
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Cb) Step function. 



j 


(c) Broken -line function. 

Fi 0 ure 5.- Typical func+ions for idealized sections 
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Figure 6.- Shear flows at jth stiffener. 
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Figure 7.- Forces acting on a differential 

length ot a stiffener. 
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(b) An+isymme+rical case. 


Figure 8.~ Modified warping s+iffness 
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(b) Antisymmetrical case - 

Figure 9. ~ Stiff ness modifications 
for a central flange. 
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(b) Double (c) Single 

(a) Stiffnesses. symmetry. symmetry 

about y-axis . 


Figure 11." Six-flange section. 
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Cb) Double (c) Single 

(a) Stiffnesses. symmetry, symmetry 

about y-axis , 


Figure IZ.- Six-flange section with two cells. 




100 


HACA TN 2529 


Q= 

4 


2 

==Q 

1 

5 

o== 

6 

7 

8 

==^ 
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Cb) Double symmetry. 


(c) Double Antisymmetry 


Figure 13. “ Eight-flange section, 
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(a) Cb) 


Figure 14 - Trapezoidal section. 
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Figure 16 . - Geometric properties associated 

with l^tb web. 
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Figure 17. - Elementarg section wHh two cells. 
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